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ABSTRACT 


In this research, a continuous analog generalized impedance converter (GIC) 4 th 
order band pass filter is investigated in detail. Various classroom software aids such as 
MATLAB, P-SPICE and MAPLE are utilized to simulate varies circuit parameter 
changes in ideal and non-ideal GIC filter, such as network sensitivity, effects of resistor 
value, capactor value and reduction of gain bandwidth product (GBWP), on the center 
frequency and Q factor. All simulated results are used to evaluate the actual circuit im¬ 
plementation prior to future GIC filter chip fabrication. 

A composite operational amplifier utilizing the BiCMOS standard operational 
amplifier that was designed and introduced in a previous thesis is presented and its im¬ 
proved performance is investigated. The composite amplifier (C20A2) is designed using 
Silvaco EXPERT and simulated with Silvaco SmartSpice. The results show the gain 
bandwidth product (GBWP), common mode rejection ratio (CMRR), and open loop gain 
are considerably improved. This sets the basic foundation for future students to incorpo¬ 
rate the newly designed composite operational amplifier into the GIC filter to further en¬ 
hance filter perfonnance. 


v 



THIS PAGE INTENTIONALLY LEFT BLANK 


vi 



TABLE OF CONTENTS 


I. INTRODUCTION.1 

A. BACKGROUND.1 

B. OBJECTIVE.1 

C. RELATED WORK.2 

D. THESIS ORGANIZATION.2 

II FILTERS.3 

A FILTER BASICS.3 

1. Filter Types.4 

a. Low-Pass Filter (LP) . 4 

b. High-Pass Filter (HP) . 4 

c. Band-Pass Filter (BP) .5 

d. Notch Filter .5 

2. Filter Quality Factor.5 

B. PASSIVE FILTERS.6 

1. How Passive Filters Work.8 

2. Mathematical Model.8 

C. ACTIVE FILTER.9 

1. Operational Amplifier.9 

a. dc Behavior . 10 

b. ac Behavior . 11 

c. Applications . 11 

d. Op-Amp Limitation . 13 

e. Internal Circuitry . 14 

f. Common Configuration . 14 

2. Type of Filters.18 

a. Chebyshev Filter . 18 

b. Butterworth Filter . 18 

c. Bessel Filter . 18 

d. Gaussian Filter . 19 

e. Sallen Key Filter. . 19 

f. Elliptic Filters . 20 

3. Advantages and Disadvantages.20 

D. GIC FILTER.21 

III. VLSI DESIGN THEORY.25 

A. INTRODUCTION [10],[11].25 

B. VLSI DESIGN FLOW [10],[11].27 

C. FABRICATION PROCESS FLOW - BASIC STEPS [12].28 

D. LAYOUT DESIGN RULES [11].31 

E. LAYOUT OF CMOS INVERTER [12].32 

F. LAYOUT OF CMOS NAND AND NOR GATES [12].35 

IV. BICMOS COMPOSITE OPERATIONAL AMPLIFIER SIMULATION 

RESULTS.37 

vii 












































A. BACKGROUND.37 

B. OPEN LOOP GAIN, 3 DB FREQUENCY AND Q P FUNCTION.39 

C. COMPOSITE OPERATIONAL AMPLIFIER BANDWIDTH 

IMPROVEMENT.42 

D. COMPOSITE OPERATIONAL AMPLIFIERS SENSITIVITIVES.42 

E. BLOCK DIAGRAM OF AUTHOR’S STANDARD AND 

COMPOSITE DESIGN USED IN THESIS.43 

F. SCHEMATIC DIAGRAM OF AUTHOR’S STANDARD AND 

COMPOSITE DESIGN USED IN THESIS.44 

G. SIMULATION OF AUTHOR’S STANDARD AND COMPOSITE 

DESIGN USED IN THESIS.45 

1. Open Loop Gain.46 

2. Closed Loop Gain.48 

3. Slew Rate.49 

4 Common Mode Rejection Ratio (CMRR).50 

H. SUMMARY OF SIMULATION PERFORMANCE INDICTOR.51 

V. GIC FILTER SIMULATION RESULTS.53 

A. DESIGN DETAILS.53 

B. IDEAL AND NON IDEAL TRANSFER FUNCTIONS.55 

C. SENSITIVITIES ANALYSIS FOR S T(S ) R , S T(S ) C , S mo R AND S QP R .55 

1. St( S ) R = S(Numerator, R) - S(Denominator, R).55 

2. S T(s) c = S(Num,C) - S(Deno,C).56 

3. S roo R = S(Num,r) - S(Deno,r).57 

4. S Qp r = S(Num,R) - S(Deno,R).58 

D. PLOTTING THE IDEAL AND NON-IDEAL TRANSFER 

FUNCTION USING MATLAB AND P-SPICE SOFWARE.58 

1. Comparing the Ideal Values with MATLAB Program Values.60 

2. Comparing the MATLAB Program Values and P-SPICE 

Values.60 

E. ANALYSIS OF THE EFFECT OF 10% VARIATION ON VALUE 

OF RESISTOR ON fo AND Qp OF THE NON-IDEAL TRANSFER 
FUNCTION.61 

1. Comparing the Values in Both MATLAB Programs (with and 

without 10% Increase in Resistor Value).62 

2. Comparing the MATLAB Simulation Program Values (with 
10% Increase in Resistor Value) and P-SPICE Values (with 

10% Increase in Resistor Value).62 

F. ANALYSIS OF THE EFFECT OF 10% VARIATION ON VALUE 

OF CAPACITOR ON fo AND Qp OF THE NON-IDEAL 
TRANSFER FUNCTION.63 

1. Comparing the Values in Both MATLAB Programs (with and 

without 10% Increase in Capacitor Value).64 

2. Comparing the MATLAB Program Values (with 10% Increase 
in Capacitor Value) and P-SPICE Values (with 10% Increase 

in Capacitor Value).64 

viii 
































G. ANALYSIS OF THE EFFECT OF 10% VALUE INCREASE IN 

RESISTOR AND 50% DECREASE IN GBWP OF THE NON¬ 
IDEAL TRANSFER FUNCTION.64 

1. Comparing the Values in Both MATLAB Programs (with 10% 

Increase in Resistor Value for 100% and 50% Reduction of 
GBWP).66 

2. Comparing the MATLAB Program Values (with 10% Increase 
in Resistor Value for 50% Reduction of GBWP) and P-SPICE 
Values (with 10% Increase of Resistor Value for 50% 

Reduction of GBWP).66 

H. ANALYSIS OF THE EFFECT OF 10% VALUE INCREASE IN 

CAPACITOR AND 50% DECREASE IN GBWP OF THE NON¬ 
IDEAL TRANSFER FUNCTION.66 

1. Comparing the Values in Both MATLAB Programs (with 10% 

Increase in Capacitor Value for 100% and 50% Reduction of 
GBWP).68 

2. Comparing the MATLAB Program Values (with 10% Increase 

in Capacitor Value for 50% Reduction of GBWP) and P- 
SPICE Values (with 10% Increase in Capacitor Value for 50% 
Reduction of GBWP).68 

I. SUMMARY OF SIMULATION RESULTS.68 

VI. CONCLUSIONS AND RECOMMENDATIONS.73 

APPENDIX A. OPEN LOOP.75 

A. OPEN LOOP - LEE’S CONFIGURATION [6], SPICE NETLIST.75 

B. OPEN LOOP - C20A2 CONFIGURATION SPICE NETLIST.76 

APPENDIX B. CLOSED LOOP.79 

A. CLOSED LOOP - LEE’S CONFIGURATION [6], SPICE NETLIST ..79 

B. CLOSED LOOP - C20A2 CONFIGURATION SPICE NETLIST.80 

APPDENDIX C. SLEW RATE.83 

A. SLEW RATE - LEE’S CONFIGURATION [6], SPICE NETLIST.83 

B. SLEW RATE - C20A2 CONFIGURATION SPICE NETLIST.84 

APPENIDX D. CMRR.87 

A. CMRR - LEE’S CONFIGURATION [6], SPICE NETLIST.87 

B. CMRR - C20A2 CONFIGURATION SPICE NETLIST.88 

APPENDIX E. DERIVATION OF IDEAL TRANSFER FUNCTION FOR 

4TH-ORDER GIC BANDPASS FILTER.91 

APPENDIX F. IDEAL GIC FILTER MATLAB PROGRAM.95 

APPENDIX G. NON- IDEAL GIC FILTER MATLAB PROGRAM.97 

APPENDIX H. NON-IDEAL MATLAB PROGRAM (NON IDEAL CASE 

WITH 10% INCREASE IN R).99 

APPENDIX I. NON-IDEAL MATLAB PROGRAM (NON IDEAL CASE 

WITH 10% INCREASE IN C).103 


IX 



























APPENDIX J. NON-IDEAL MATLAB PROGRAM (NON-IDEAL CASE 

WITH 10% INCREASE IN R AND 50% DECREASE OF GBWP).107 

APPENDIX K. NON-IDEAL MATLAB PROGRAM (NON-IDEAL CASE 

WITH 10% INCREASE IN C AND 50% DECREASE OF GBWP).Ill 

LIST OF REFERENCES.115 

INITIAL DISTRIBUTION LIST.117 


x 







LIST OF FIGURES 


Figure 1. Quality Factor for a Band-Pass Filter (From [3]).6 

Figure 2. Op-Amp Symbol.10 

Figure 3. Voltage Divider.12 

Figure 4. Inverting Amplifier.15 

Figure 5. Non-Inverting Amplifier.15 

Figure 6. Voltage Follower.16 

Figure 7. Summing Amplifier.17 

Figure 8. Intertegrator.17 

Figure 9. Differentiator.17 

Figure 10. Basic GIC Filter (From [2]).21 

Figure 11. Prominent Trends in Information Service Technologies. (From [10]).25 

Figure 12. Evolution of Integration Density and Minimum Feature Size, as Seen in 

the Early 1980’s (From [10]).26 

Figure 13. VLSI Design Flow (From [10]).27 

Figure 14. Process Steps Required for Patterning of Silicon Dioxide (From [12]).29 

Figure 15. The Result of a Single Lithographic Patterning Sequence on Silicon 

Dioxide, without Showing the Intennediate Steps. (From [12]).31 

Figure 16. Design Rule Constraints which Determine the Dimensions of a Minimum- 

Size Transistor (From [12]).33 

Figure 17. Placement of One nMOS and One pMOS Transistor (From [12]).34 

Figure 18. Complete Mask Layout of the CMOS Inverter (From [12]).34 

Figure 19. Layouts of a CMOS NOR2 Gate and a CMOS NAND2 Gate. (From [10]) ...35 

Figure 20. C20A1 (From [15]).38 

Figure 21. C20A2 (From [15]).38 

Figure 22. C20A3 (From [15]).39 

Figure 23. C20A4 (From [15]).39 

Figure 24. Block Diagram of Standard Op-Amp.43 

Figure 25. Block Diagram of Composite Op-Amp C20A2.44 

Figure 26. Circuit Layout of Lee Configuration.44 

Figure 27. Symbolic Diagram of C20A2.45 

Figure 28. Schematic Diagram of C20A2 Design.45 

Figure 29. Open Loop Transfer Function of Lee’s Configuration B.46 

Figure 30. Open Loop Transfer Function of C20A2.47 

Figure 31. Open Loop Transfer Function - Comparison.47 

Figure 32. Closed Loop Gain Configuration.48 

Figure 33. Closed Loop Gain of Lee’s Configuration.49 

Figure 34. Closed Loop Gain of C20A2.49 

Figure 35. Slew Rate of Lee’s Configuration.50 

Figure 36. Slew Rate of C20A2.50 

Figure 37. CMRR of Lee’s Configuration and C20A2.51 

Figure 38. Circuit Schematic.54 

Figure 39. Circuit Layout.54 


xi 









































Figure 40. MATLAB Ideal Transfer Function.59 

Figure 41. MATLAB Non-Ideal Transfer Function.59 

Figure 42. P-SPICE Frequency Response Plot.60 

Figure 43. MATLAB Simulation (Non-Ideal Case with 10% Increase in Resistor 

Value).61 

Figure 44. P-SPICE Frequency Response Plot (with 10% Increase in Resistor Value)...62 
Figure 45. MATLAB Simulation (Non-Ideal Case with 10% Increase in Capacitor 

Value).63 

Figure 46. P-SPICE Frequency Response Plot (with 10% Increase in Capacitor 

Value).64 

Figure 47. MATLAB Plot (Non-Ideal Case with 10% Increase in Resistor Value and 

50% Decrease of GBWP).65 

Figure 48. P-SPICE Frequency Response Plot (with 10% Increase in Resistor Value 

and 50% Decrease of GBWP).66 

Figure 49. MATLAB Plot (Non-Ideal Case with 10% Increase in Capacitor Value 

and 50% Decrease of GBWP).67 

Figure 50. P-SPICE Frequency Response Plot (with 10% Increase in Capacitor 

Value and 50% Decrease of GBWP).68 

Figure 51. MATLAB Plot: Center Frequency = 49.338 kHz, Qp = 10.33.70 

Figure 52. P-SPICE Frequency Response: Center Frequency = 47.853 kHz, Qp = 

10.655.71 

Figure 53. Actual Frequency Response: Center Frequency = 50.560 kHz, Qp = 9.711 ..71 


xii 















LIST OF TABLES 


Table 1. Ideal Frequency Response of the Four Basic Filter Types. (From [2]).3 

Table 2. Non-Ideal Filters Response [From 3].7 

Table 3. Transfer Function of Four Basic Filter Types. (From [3]).9 

Table 4. Admittance Selection and Transfer Functions (From [1]).22 

Table 5. C20As Open Loop Gain Input-Output Relationships (From [15]).40 

Table 6. C20A 3dB Frequency and Qp Functions (From [15]).41 

Table 7. Stability Criteria for C20As (From [15]).41 

Table 8. Summary of Perfonnance Indictor.51 

Table 9. Summary of Component Values.69 

Table 10. Summary of Calculated Values.70 


xiii 













THIS PAGE INTENTIONALLY LEFT BLANK 


xiv 



ACKNOWLEDGMENTS 


The author would first like to thank Professor Michael for his enthusiasm, guid¬ 
ance and support throughout the duration of this project. He also wants to thank Professor 
Fouts for always answering the persistent questions. The author’s wife deserves special 
recognition for all her support during the past one and a half years. Without her help and 
understanding regarding the late hours at work, this would have been a much more diffi¬ 
cult ordeal. 


xv 



THIS PAGE INTENTIONALLY LEFT BLANK 


xvi 



EXECUTIVE SUMMARY 


Electronic filters and their circuit applications can be found in almost any vital 
military systems hardware. They play an essential role in all comminaction systems, avia¬ 
tion electronis, GPS receivers, as well as in most digital systems. The use of an analog 
filter can eliminate some of the errors that occur during analog-to-digital conversion, 
digital processing, and digital-to-analog conversion. This thesis presents a continuous 
analog generalized impedance converter (GIC) 4th-order band pass Filter. MATLAB, P- 
SPICE and MAPLE software are utilized to simulate various circuit parameter changes in 
the Ideal and Non-Ideal GIC filter, such as network sensitivity, effects of resistor value, 
capactor value and reduction of the gain bandwidth product (GBWP), on the center fre¬ 
quency and Q factor. All simulated results are used to evaluate the actual circuit imple¬ 
mentation prior to future GIC filter chip fabrication. The center frequency and Q factor 
are affected by capacitor values changes. The reduction of GBWP reduces the gain, cen¬ 
ter frequency and Q factor. All the P-SPICE program values, when compared to 
MATLAB program values, are well within a 10% error. 

This thesis also researched the composite amplifier originally introduced by Pro¬ 
fessors Sherif Micheal and Wasty B. Mikhael at West Virginia University. The composite 
amplifier (C20A2) has shown that the GBWP, common mode rejection ratio (CMRR) 
and open loop gain have improved considerably. A composite operational amplifier, util¬ 
izing the BiCMOS standard operational amplifier that was designed and introduced in a 
previous thesis, is presented and its improved performance was investigated. This sets the 
basic foundation for future research to incorporate the newly designed composite opera¬ 
tional amplifier into the GIC filter to enhance filter performance further. 
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I. 


INTRODUCTION 


A. BACKGROUND 

Integrated circuit technology has seen a great deal of change in the past 30 years. 
Many advances have been made with digital circuits. Digital circuitry of the 2000’s is ca¬ 
pable of performing very complex computation functions. However, as good as digital 
circuits have become, they must still interface with the analog world. This is particularly 
true, especially in advanced military systems where the use of digital hardware is becom¬ 
ing more popular. 

Filters are the main components in many electronic circuits. Electronic filters 
were initially constructed from resistors, capacitors, and inductors. They were complex, 
costly, and sensitive to component variation. Operational amplifiers (op-amps) allowed 
the construction of active filters, which can eliminate a costly passive element, the induc¬ 
tor. 

In this thesis, the generalizedd impedance converter (GIC) 4th-order band pass 
filter presented is one of the common analog filters that can be used to aid in the interface 
between the digital and analog worlds. The GIC filter has proven itself a robust design 
that is highly insensitive to component variation. 

B. OBJECTIVE 

The primary objective of the research was to investigate the continuous analog 
Generalized Impedance Converter 4th-order band pass filter, ideal and non-ideal cases, 
in detail with software aids such as MATLAB, P-SPICE and MAPLE. Simulated results 
were used to compare the actual circuit implementation prior to future GIC filter chip 
fabrication. 

The secondary objective was to investigate the previously designed BiCMOS 
standard operational amplifier, introduced in a former research [7], to be utilized to im¬ 
prove the performance of the composite operational amplifier. This sets the basic founda¬ 
tion for future students to incorporate the newly designed composite operational amplifier 
into the GIC filter to enhance filter perfonnance further. 
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C. RELATED WORK 

Previous theses have focused on the design of a standard operational amplifier. 
The latest design was developed by LCDR Paul R. Milne [7], who replaced the CMOS 
operational amplifier with a BiCOMS version. This thesis uses LCDR Paul R. Milne’s 
design as a basis to improve it as a composite operational amplifier. 

D. THESIS ORGANIZATION 

All the basic analog material are included assuming the reader has no analog filter 
background. Chapter II discusses the basic types of filters and components used. It out¬ 
lines the two main type of filters, passive and active, as well as the advantages and disad¬ 
vantages of each type. The ideal and non-ideal characteristics, transfer function, and fil¬ 
ter types of the GIC filter are also discussed. Chapter III discusses the VLSI design flow, 
the fabrication process, and some basic layouts for the common gates. 

Chapter IV introduces the composite amplifier and investigates the author’s simu¬ 
lated circuits and compares them to the standard operational amplifier. Chapter V focuses 
on a continuous analog Generalizedd Impedance Converter (GIC) 4 th order band pass fil¬ 
ter. The effect of varies parameter such as network sensitivity, effects of resistor values, 
capactor values and GBWP on center frequency and Q factor are presented. All simulated 
results are used to evaluate the actual circuit implementation prior to future GIC filter 
chip fabrication. Chapter VI summarizes the results and makes further recommendations 
for future research topics. 
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II FILTERS 


A FILTER BASICS 

In this chapter a brief background on filters is presented. A filter is a device that 
consists of a group of components including resistors, capacitors, inductors and some¬ 
times active devices such as operational amplifiers. The primary purpose of most filters is 
to modify the amplitude response of a circuit to a signal so that certain frequencies are at¬ 
tenuated or blocked, while others pass through unchanged. [1] 

Most filters are placed in a group based on frequency ranges that are attenuated or 
passed. Generally, the ideal transmission characteristics of the four major filters are the 
low-pass (LP) filter (which passes low frequencies), the high-pass (HP) filter (which 
passes high frequencies), the band-pass (BP) filter that passes a limited range of midband 
frequencies, and the notch (N) filter (which is a type of band-stop filter that acts on a par¬ 
ticularly narrow range of midband frequencies). Table 1 shows these characteristics. 
These idealized characteristics are known as brick-wall type responses because of their 
vertical edges. [2] 



Table 1. Ideal Frequency Response of the Four Basic Filter Types. (From [2]) 
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1 . 


Filter Types 

a. Low-Pass Filter (LP) 

A low-pass filter is a filter that allows low frequencies to pass through but 
attenuates high frequencies [2 &3]. Low-pass filters are used to block unwanted high- 
frequency signals, while passing the lower frequencies. The low frequencies to be filtered 
out are relative to the unwanted higher frequencies and, therefore, do not have a defini¬ 
tive range. The frequencies that are cut vary from filter to filter. They are also called 
treble cut filters. They are the opposite of high-pass filters. 

One example of low-pass filters is a physical barrier acting as a low-pass 
filter for waves. When music is playing in another room, the low notes are easily heard, 
while the high notes are largely filtered out. Similarly, very loud music played in one car 
is heard as a low throbbing by occupants of other cars since the closed vehicles and air 
gap function as a very low frequency low-pass filter. 

Low-pass filters are also used in subwoofers and other types of loud¬ 
speakers to block high pitches that they cannot efficiently reproduce. Radio transmitters 
use low-pass filters to block hannonic emissions, which might cause interference with 
other communications. 

b. High-Pass Filter (HP) 

A high-pass filter passes ’high’ frequencies fairly well, but attenuates ’low' 
frequencies. Hence, it is useful as a filter to block any unwanted low frequency compo¬ 
nents of a complex signal while passing the higher frequencies, which is better known as 
a bass-cut filter [2 &3]. Of course, the meanings of'low' and 'high' frequencies are rela¬ 
tive. 

The simplest high-pass filter consists of a capacitor in series with the sig¬ 
nal path in conjunction with a resistor parallel to the signal path. The resistance times the 
capacitance (RC) is the time constant and its reciprocal is the cutoff frequency, at which 
the output voltage is 70.7% of the input. Such a circuit might be used in combination with 
a tweeter and a loudspeaker. 
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c. Band-Pass Filter (BP) 

A band-pass filter passes a limited range of frequencies, and possibly is 
created as a combination of a low-pass filter and a high-pass filter [2 &3]. 

A band-pass filter will allow only frequencies to pass that are in a previ¬ 
ously specified range. For example, in atmospheric sciences, it is common to perform 
band-pass filtering of the data with a period range of, say 3 to 10 days, so that only 
cyclones remain as fluctuations in the fields. 

A band-pass filter is an electronic circuit that passes through signals be¬ 
tween two given frequencies. For example, an ideal band-pass filter might allow all sig¬ 
nals to pass above 30 Hz but below 100 Hz. All signals outside this range are attenuated. 

In practice, most band-pass filters do not attenuate frequencies just out¬ 
side the desired frequency range completely. 

d. Notch Filter 

A notch filter, also called a narrow band-stop filter or T-notch filter, is an 
electronic filter typically used when the high frequency and the low frequency are less 
than 1 to 2 decades apart. In other words, the high frequency is less than 10 to 20 times 
the low frequency [2 &3]. 

In the audio band, a notch filter uses high and low frequencies that may be 
only semitones apart [2]. 

2. Filter Quality Factor 

The quality factor, or “Q”, is a parameter used to describe the selectivity perform¬ 
ance of a filter. For a first-order filter, the Q parameter relates the distance of the filter 
pole to the ja>— axis. For a more selective filter, the Q value must be high. A high Q im¬ 
plies that the poles of the filter have to be closer to the jco- axis. The Q of a filter can be 
calculated by taking the inverse of the normalized bandwidth of the filter. Calculations 
for normalization can be achieved by using the half power (3dB) point of the filter. Fig¬ 
ure 1 illustrates the Q calculation. [3] 

Q = ~^~ ( 2 - 1 ) 

(o 1 —co 2 
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COq CO 

co x co 2 

Figure 1. Quality Factor for a Band-Pass Filter (From [3]) 

B. PASSIVE FILTERS 

The simplest electronic filters are based on combinations of resistors, inductors 
and capacitors. Since resistance has the symbol R, inductance the symbol L and 
capacitance the symbol C, these filters exist in so-called RC, RL, LC and LCR varieties. 
All these types are collectively known as passive filters, because they are activated by the 
power in the signal and not by an external power supply. 

The figures in Table 2 show the RLC network and the non-ideal response charac¬ 
teristics of the four basic types of filters. Although ideal filters cannot be realized, their 
concepts aid in filter analysis. Non-ideal filters, unlike their ideal counterparts, do not 
transition perfectly but have a transition band from the pass-band to the stop-band. 
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Table 2. Non-Ideal Filters Response [From 3] 
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1. How Passive Filters Work 

Inductors block high-frequency signals and conduct low-frequency signals, while 
capacitors do the reverse. A filter in which the signal passes through an inductor, or in 
which a capacitor provides a path to ground, therefore transmits low-frequency signals 
more strongly than high-frequency signals and is a low-pass filter. If the signal passes 
through a capacitor, or has a path to ground through an inductor, then the filter transmits 
high-frequency signals more strongly than low-frequency signals and is a high-pass fil¬ 
ter. Resistors on their own have no frequency-selective properties, but are added to in¬ 
ductors and capacitors to determine the time-constants of the circuit, and therefore, the 
frequencies to which it responds [2 &3]. 

At very high frequencies above approximately 100 megaHertz, sometimes the in¬ 
ductors consist of single loops or strips of sheet metal, and the capacitors consist of adja¬ 
cent strips of metal. These are called stubs. It is possible to add other components to LC 
filters to make them more precise. 

2. Mathematical Model 

Filter characteristics are specified using a mathematical model, the transfer func¬ 
tion. The exponential degree of the denominator is known as the order of the filter. To 
achieve all the basic filter types, the transfer function must be at least second order. 

Based on the mathematical model, the actual physical filters using the resistors 
(. R ), capacitor (C) and inductor (L) can be realized. Table 3 shows the generalizedd trans¬ 
fer functions for each filter type. The center frequency is given by co p , while Q gives 

the quality factor, and a j is a constant. [3] 
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Filter Type 

Transfer Function 

Low-Pass 

T(s)= - & -- 

Qp 

— 

High-Pass 

t(s)= — sd- 

s ^+^JL s + co * 

Qr " 

Band-Pass 

T (s) = - ^ - 

Qp 

Notch 

1 r(s)= ^ 

2 eo , 

s +—-S + 0)„ 

Qr P 


Table 3. Transfer Function of Four Basic Filter Types. (From [3]) 


C. ACTIVE FILTER 

Active filters are a special kind of electronic filter. They use additional energy 
than that of the signal to filter the signal. They are commonly constructed using an 
operational amplifier. 

The topologies correspond to R and C components along a retro alimented loop 
around the operational amplifier. These can have high Q, and achieve resonance without 
the use of inductors. However, their upper frequency limit is lower than that of a passive 
filter. Some common kinds of filters are Chebyshev, Butterworth, Bessel, Gaussian, 
Sallen and Key, and Elliptic. 

1. Operational Amplifier 

An operational amplifier or op-amp is an electronic circuit module normally built 
as an integrated circuit, which has a non-inverting input (+), an inverting input (-) and 
one output. The output voltage is the difference between the + and - inputs multiplied by 
the open-loop gain V out = (V+ - V-) * G op enioop. Since op-amps have uniform parameters 
and often standardized packaging as well as standard power supply needs, they help in 
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designing an application quickly. Figure 2 illustrates a typical circuit symbol for an op- 
amp, where V+is the non-inverting input, V- the inverting input, V out the output, Vs+ the 
positive power supply and Vs- the negative power supply. 

v R+ 



Figure 2. Op-Amp Symbol 


The power supply pins (Vs+ and Vs-) can be labeled in many different ways. For 
field-effect-transistor (FET) based op-amps, the positive, common drain supply is la¬ 
beled Vdd and the negative, common source supply is labeled Vss- For a bipolar junction 
transistor (BJT) based opamps, the Vs+ pin becomes Vcc and Vs- becomes V EE . They are 
also sometimes labeled Vcc+ and Vcc-, or even V+ and V-, in which case the inputs 
would be labeled differently. The function remains the same. Often these pins are not in¬ 
cluded in the diagram for purposes of clarity, and the power configuration described or 
assumed from the circuit. 

The input pin polarity is often reversed in diagrams for clarity. In this case, the 
power supply pins remain in the same position. The more positive power pin is always on 
the top, and the more negative on the bottom. The entire symbol is not flipped; only the 
inputs. 

a. dc Behavior 

Open-loop gain is defined as the amplification from input to output with¬ 
out any feedback applied. For most practical calculations, the open-loop gain is assumed 
to be infinite. This allows the gain in the application to be set simply and exactly by using 
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negative feedback. Of course, theory and practice differ, since op-amps have limits that 
the designer must keep in mind and sometimes adjust the work. 

b. ac Behavior 

The op-amp gain calculated at dc does not apply at higher frequencies. 
This effect is due to limitations within the op-amp itself, such as its finite bandwidth, and 
to the ac characteristics of the circuit in which it is placed. The best-known impediment 
in designing with op-amps is the tendency for the device to resonate at high frequencies, 
where negative feedback changes to positive feedback due to internal parasitics capaci¬ 
tances. 

c. Applications 

The op-amp can be used for various applications such as audio and video 
pre-amplifiers, voltage comparators, differential amplifiers, differentiators and integra¬ 
tors, filters, precision rectifiers, voltage and current regulators, analog calculators, ana- 
log-to-digital converters, digital-to-analog converters and many more. 

The operational amplifier is thus called because it performs mathematical 
operations by using voltage as an analog of another quantity. This is the basis for the ana¬ 
log computer. As mentioned before, the generic op-amp differential gain is given by: 

Vout = G (V+ - V-). (2.2) 

G is the open-loop gain of the op-amp. The inputs are assumed to have 
very high impedance, and negligible current will flow into or out of the inputs. Op-amp 
outputs have very low source impedance. If the output is connected to the inverting input, 
after being scaled by a voltage divider, K = Ri / (Ri + R 2 ) as shown in Figure 3. 
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out 

Figure 3. Voltage Divider 

Vout = G(V m -KV out ). (2.3) 

The result is a linear amplifier with gain: 

Vout / Vin = G / (1 + G K). (2.4) 

If G is very large, this comes close to 1 / K = 1 + (R 2 / Ri). 

This negative feedback connection is the most typical use of an op-amp, 
but many different configurations are possible, making it one of the most versatile of all 
electronic building blocks. 

When connected in a negative feedback configuration, the op-amp will 
tend to produce whatever voltage is necessary to make the input voltages equal. This, and 
the high input impedance, are sometimes called the two “golden rules” of op-amp design 
for circuits that use feedback. No current will flow into the inputs and the input voltages 
will be equal to each other. The exception is if the voltage required is greater than the op- 
amp's supply, in which case, the output signal stops at a value very close to the voltages, 
V s + or V s - 
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Most integrated single, dual and quad op-amps circuits available have a 
standardized pin out, which permits one type to be substituted for another without wiring 
changes. A specific op-amp may be chosen for its open loop gain, bandwidth, noise per¬ 
formance, input impedance, power consumption, or a compromise between any of these 
factors. 

Historically, the first integrated op-amp to become widely available was 
the Fairchild UA-709, in the late 1960’s, but this was rapidly superseded by the much 
better perfonning 741, which is easier to use, and probably ubiquitous in electronics. 

Most manufacturers produce a version of this classic chip. The 741 is a bipolar design, 
and by modern standards has fairly average perfonnance. Better designs based on the 
FET arrived in the late 1970’s, and MOSFET versions in the early 1980’s. It is possible 
to substitute many of these more modern devices into an older 741-based circuit and 
work with no other changes to provide better performance. 

d. Op-Amp Limitation 

Although the design of most op-amp circuits relies on the “golden rules” 
above, designers should also be aware that no real op-amp can match these characteris¬ 
tics exactly. Listed below are some of the limitations of real op-amps, as well as how this 
affects circuit design. 

(1) DC Imperfection. There are five main DC imperfections as 
discussed in this section. Finite gain is the effect that is most pronounced when the over¬ 
all design attempts to achieve gain close to open loop gain of the op-amp. Finite input 
resistance puts an upper bound on the resistances in the feedback circuit. Nonzero output 
resistance is important for low resistance loads. Except for very small voltage outputs, 
power considerations usually come into play first. Input bias current, a small amount of 
current, typically ~10nA, into the input pins, is required for proper operation. This effect 
is aggravated because this current is mismatched slightly between the input pins (i.e., in¬ 
put offset current). This effect is usually important only for very low power circuits. Input 
offset voltage is the op-amp that will produce an output even when the input pins are at 
exactly the same voltage. For circuits which require precise DC operation, it is necessary 
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to compensate for this effect. Most commercial op-amps provide an offset pin for this 
purpose. 

(2) AC Imperfections. Two AC imperfections are discussed. 
The first is a finite bandwidth which all amplifiers possess. However, this is more pro¬ 
nounced in op-amps which use frequency compensation to avoid unintentionally produc¬ 
ing positive feedback. The second is the input capacitance, which is the most important 
for high frequency operation. 

(3) Nonlinear Imperfections. Two nonlinear imperfection are 
discussed. Saturation refers to output voltage that is limited to a peak value slightly less 
than the power supply voltage. The slew rate is the limit of the rate of change of the out¬ 
put voltage. 

(4) Power Considerations. Most op-amps are not designed for 
high power operation. For high power op-amp circuits, an op-amp specifically designed 
for that purpose must be used. Short circuit protection is more a feature than a limitation, 
although it does limit the design. Most commercial op-amps shut off when the load resis¬ 
tance is below a specified level. 

e. Internal Circuitry 

Although it is useful and easy to treat the op-amp as a black box with a 
perfect input/output characteristic, it is important to understand the inner workings, so 
that it is possible to deal with problems arising due to internal parasitic capacitances. 

Although designs vary between products and manufacturers, all op-amps 
have basically the same internal structure, which consists of three stages. The differential 
amplifier input stage provides low noise amplification, high input impedance and drives a 
current mirror load. The voltage amplifier provides high gain. The output amplifier out¬ 
put stage provides high current driving capability, low output impedance, current limiting 
and short circuit protection. 

/. Common Configuration 

The resistors used in these configurations are typically in the k Q range. 
The <1 kQ range resistors cause excessive current flow and possible damage to the de¬ 
vice. The >1 Mfi range resistors cause excessive thermal noise and undesired bias cur- 
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rents. The Z out for all amplifiers is ideally 0 Q . Realistically, it is 1 Q to 1 kQ , depend¬ 
ing on the device. 

(1) Inverting Amplifier. Figure 4 shows the configuration for 
an inverting amplifier. It inverts and amplifies a voltage, which multiplies by a negative 
constant, and because V- is a virtual ground. 

V out = V in (Rf / Rin) (2.5) 

Z in = R in . (2.6) 


R f 



Figure 4. Inverting Amplifier 

(2) Non-Inverting Amplifier. Figure 5 shows the configura¬ 
tion for a non-inverting amplifier. It amplifies a voltage, which multipies by a constant 

greater than 0 but less than 1. Realistically, the input impedance of the circuit is the input 
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impedance of the op-amp itself, 1 MQ to 10 Q. 


Vout = V.n(l+R2/Rl) (2.7) 

Zin = °o (2.8) 



Figure 5. Non-Inverting Amplifier 
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(3) Voltage Follower. Figure 6 shows the configuration for a 
voltage follower, which is used as a buffer to eliminate loading effects or to interface im¬ 
pedances connecting a device with a high source impedance to a device with a low input 
impedance. 

Vout = Vin (2.9) 

Zj n = oo (2.10) 



Figure 6. Voltage Follower 

(4) Summing Amplifier. Figure 7 shows the configuration for 
a summing amplifier. It sums several (weighted) voltages and the output is inverted. The 
input impedance is Z n = R n for each input and V- is a virtual ground). For independent 
Ri, R?, ... R n , 

V = - Rf (Vi / Ri + V 2 / R 2 + ... + Vn / Rn) (2.11) 

For Ri = R 2 = ... = R n 

V = -(Rf/R 1 )(V 1 +V 2 + ... + V„) (2.12) 

For Rf = Ri = R 2 = ... = R n 

V = -(V!+V 2 +...+Vn). (2.13) 


16 




Figure 7. Summing Amplifier 

(5) Integrator and Differentiator. Figures 7 and 8 show the 
configuration for an integrator and differentiator, respectively. The integrator integrates 
the signal over time. The differentiator differentiates the signal over time. 


c 



Figure 8. Intertegrator 

R 



Figure 9. Differentiator 
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2. Type of Filters 

Several common types of filters are Chebyshev, Butterworth, Bessel, Gaussian, 
Sallen and Key, and Elliptic. 

a. Chebyshev Filter 

The Chebyshev Type I filter is the filter type that results in the sharpest 
pass band cut off and contains the largest group delay. The most notable feature of this 
filter is the ripple in the pass band magnitude. 

The pass band attenuation of the standard Chebyshev Type I filter is de¬ 
fined to be the same value as the pass band ripple amplitude. However, Filter Solutions 
allows the user the option of selecting any pass band attenuation that will define the fil¬ 
ter’s cut off frequency. 

The Chebyshev Type II filter, also known as the Inverse Chebyshev filter, 
contains a Butterworth style, or maximally flat, pass band, a moderate group delay, and 
an equiripple stop band. Like the Butterworth filter, the pass band attenuation of the Che¬ 
byshev Type II filter is defined to be -3.01 dB. However, Filter Solutions allows the user 
the option of selecting any pass band attenuation in dB’s that will define the filter’s cut 
off frequency. 

b. Butterworth Filter 

The Butterworth filter is the filter type that results in the flattest pass band 
and contains a moderate group delay. A standard Butterworth filter's pass band attenua¬ 
tion is -3.01dB. However, Filter Solutions allows the user the option of selecting any 
pass band attenuation in dB's that will define the filter’s cut off frequency. 

Filter Solutions also offers the user the option of placing user-defined ze¬ 
ros in the stop band. Such a filter with stop band zeros is no longer a true Butterworth Fil¬ 
ter, but is still in the Maximally Flat filter family. 

c. Bessel Filter 

The Bessel filter's distinguishing characteristic is the near constant group 
delay throughout the pass band of the low pass filter. Filter Solutions normalizes the Bes¬ 
sel filter such that the prototype high frequency attenuation matches the Butterworth fil¬ 
ter. The pass band attenuation of the Bessel filter increases with the order of the filter 
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when this normalization is applied. However, Filter Solutions allows the user the option 
of selecting the desired pass band attenuation in dB’s. 3dB attenuation is a popular choice 
for some. 

d. Gaussian Filter 

The Gaussian filter is the filter type that results in the most gradual pass 
band roll-off and the lowest group delay. The step response of the Gaussian filter 
NEVER overshoots the steady state value. As the name states, the Gaussian filter is de¬ 
rived from the same basic equations used to derive the Gaussian distribution. The signifi¬ 
cant characteristic of the Gaussian filter is that the step response contains no overshoot at 
all. 

Filter Solutions normalizes the Gaussian filter such that the prototype high 
frequency attenuation matches the Butterworth filter. The pass band attenuation of the 
Gaussian filter increases with the order of the filter when this nonnalization is applied. 
However, Filter Solutions allows the user the option of selecting the desired pass band at¬ 
tenuation. 3dB attenuation is a popular choice. 

e. Sallen Key Filter 

A Sallen Key filter is a two poles filter topology. It is sometimes called 
VCVS, which stands for voltage controlled voltage source. It is available in low pass, 
high pass, bandpass, and notch versions, although the notch variety owes more to the 
Twin T topology. The bandpass version of the Sallen Key topology is not recommended, 
because the input resistor tends to be a very low value. The Sallen Key topology is often 
seen in a unity gain version, although the addition of gain resistors is an obvious en¬ 
hancement. The designer is cautioned, however, that changing the gain also changes 
other filter characteristics. As in the other single opamp topologies, changing gain will 
also affect the frequency and filter type (Butterworth, Chebyshev, and Bessel). 

The Sallen Key topology is suited for operation from a single supply. Cur¬ 
rent feedback amplifiers can be used, with the restriction that the connection from the op 
amp output to inverting input must be the recommended feedback resistor value. This is 
true whether or not this resistor is part of a gain stage. It is not possible to use the Sallen 
Key topology with fully differential amplifiers. 
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/. Elliptic Filters 

The Elliptic filter contains a Chebyshev Type I style equiripple pass band, 
an equipped stop band, a sharp cut off, high group delay, and the greatest stop band at¬ 
tenuation. 

Like the Chebyshev Type I filter, the Elliptic pass band attenuation is de¬ 
fined to be the same value as the pass band ripple amplitude. However, Filter Solutions 
allows the user the option of selecting any pass band attenuation that will define the fil¬ 
ter’s cut off frequency. 

3. Advantages and Disadvantages 

A filter is a circuit that allows certain frequencies to pass and blocks others. This 
selective nature can be done two ways, either with passive components or with active 
components. Passive filters are completely comprised of passive elements; namely resis¬ 
tors, capacitors and/or inductors. Active filters use active devices, i.e., an op-amp, to fil¬ 
ter out unwanted signals. 

Passive filters have some important disadvantages in certain applications, how¬ 
ever. Since they use no active elements, they cannot provide signal gain. Input imped¬ 
ances can be lower than desirable, and output impedances can be higher than optimal for 
some applications, and therefore, buffer amplifiers may be needed. Inductors are neces¬ 
sary for the synthesis of the most useful passive filter characteristics, and these can be 
prohibitively expensive if high accuracy, 1% or 2% for example, small physical size, or a 
large value are required. Standard values of inductors are not very closely spaced, and it 
is diffcult to find an off-the-shelf unit within 10% of any arbitrary value. Thus, adjust¬ 
able inductors are often used. Tuning these to the required values is time-consuming and 
expensive when producing large quantities of filters. Futhermore, complex passive filters 
higher than 2nd-order can be difficult and time-consuming to design. 

On the other hand, active filters have high gain and it is easy to adjust and tune 
the frequency. There are no inductors, which will reduce the cost and size of the circuit. It 
has no loading effects. Some disadvantages of active filters include the bandwidth limita¬ 
tions, fabrication tolerances, and that it can only respond to a specific range of signal 
magnitudes. 
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D. GIC FILTER 


The generalizedd impedance converter (GIC) is highly insensitive to component 
variation [1],[4]. The GIC filter design was introduced by Mikhael and Bhattachararyya 
[2] and proved to be very insensitive to non-ideal component characteristics and varia¬ 
tions in component values [5]. Figure 10 shows the general topology of the GIC filter. 

GIC Biquads are two op-amp biquads with good high frequency perfonnance. All 
but the even notch stages are tunable. The high pass, low pass and band pass stages are 
gain adjustable. The notch and all pass stages have a fixed gain of unity. All GIC stages 
have equal capacitor values, unless a capacitor is required to adjust the gain. Notch stages 
do not rely on element value subtractions for notch quality and are thus immune from de¬ 
gradations in notch quality due to element value errors. 



Figure 10. Basic GIC Filter (From [2]) 
The transfer functions for the GIC filter are as follows [2], 

T z Yl = 

1 V t F 1 F 4 (F 5 +F 6 ) + F 2 F 3 (F 7 +F 8 ) 


T = Yi TO+W, +WrTO 

2 Vi T,F 4 (F 5 + F 6 ) + F 2 F 3 (F 7 + F s ) 
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(2.15) 







Choosing the appropriate elements and substituting the proper admittance for each 
element into the transfer functions in Equations 2.14 and 2.15 yields transfer functions 
representing the four filter types. Table 4 depicts the transfer functions and admittance se¬ 
lection for each of the eight elements. 
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Table 4. Admittance Selection and Transfer Functions (From [1]) 

These transfer functions are based on the use of ideal op-amps. Non-ideal trans¬ 
fer functions are derived and shown below [1]. 

The non-ideal transfer function for T1 is: 
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The non-ideal transfer function for T2 is: 
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.(2.17) 


where a> t and oo t are the gain bandwidth product (GBWP) of the two op amps used. 
Note that {co t —► oc), assuming ideal op-amps, the non-ideal transfer function equations 
simplify the ideal equations from Table 4. 
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III. VLSI DESIGN THEORY 


Over the last two decades, the electronics industry has experienced phenomenal 
growth, mainly due to the rapid advances in integration technologies and large-scale sys¬ 
tems design. The number of applications of integrated circuits in high-performance com¬ 
puting, telecommunications, and consumer electronics has been rising steadily at a very 
fast pace. The fast development of this field is driven by the requirement of the computa¬ 
tional power of these applications. 

A. INTRODUCTION [10],[11] 

In the area of low-bit-rate video and cellular communications, the current lead¬ 
ing-edge technologies has already provide the end-users a certain amount of processing 
power and portability. With this trend expected to continue, it will have very important 
implications for VLSI and systems design. Figure 11 gives an overview of the prominent 
trends in information technologies over the next few decades. 



Figure 11. Prominent Trends in Information Service Technologies. (From [10]) 
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The need to integrate these functions into a small system or package is increasing 
as more and more complex functions are required in various data processing and tele¬ 
communications devices. For the past three decades, the level of integration as measured 
by the number of logic gates in a monolithic chip has been steadily rising, mainly due to 
the rapid progress in processing technology and interconnects technology. The logic 
complexity per chip has been increasing exponentially over the past few years. 

Figure 12 shows the history and forecast of chip complexity, and minimum fea¬ 
ture size over time, as seen in the early 1980’s. At that time, it was expected that it would 
be around 2000 before a minimum feature size of 0.3 microns could be realized. The ac¬ 
tual development of the technology has far exceeded these expectations. A minimum size 
of 0.25 microns was readily achievable by 1995 and the integration density has also ex¬ 
ceeded previous expectations. By the end of 1994, the first 64 Mbit DRAM, and the 
INTEL Pentium microprocessor chip containing more than three million transistors were 
available, thus pushing the envelope of integration density. 



Figure 12. Evolution of Integration Density and Minimum Feature Size, as Seen in 

the Early 1980’s (From [10]) 
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B. VLSI DESIGN FLOW [10],[11] 

The VLSI design process starts with a given set of requirements. The initial de¬ 
sign is developed and tested against the requirements and improvements made when re¬ 
quirements are not met. If such improvement is either not possible or too costly, it is then 
necessary to consider the revision of requirements and its impact analysis. 

Figure 13 shows a simplified view of the VLSI design flow, taking into account 
the various representations, or abstractions of design-behavioral, logic, circuit and mask 
layout. During this process, the verification of design plays a very important role in 
every step. The failure to verify a design properly in its early phases typically causes sig¬ 
nificant and expensive re-design at a later stage, which ultimately increases the time-to- 
market. 


VLSI 
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Behavioral 

Representation 
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(Gate-level) 

Representation 


Circuit 

Representation 


Layout 

Representation 



Fabrication & Testing 


Figure 13. VLSI Design Flow (From [10]) 
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In reality, there is no truly unidirectional top-down design flow, although top- 
down design flow provides an excellent design process control. Both top-down and bot¬ 
tom-up approaches have to be combined. For example, if a chip designer defined archi¬ 
tecture without close estimation of the corresponding chip area, then it is very likely that 
the resulting chip layout exceeds the area limit of the available technology. When this 
happened, in order to fit the architecture into the allowable chip area, some functions may 
have to be removed and the design process must be repeated. This change may require 
significant modification of the original requirements. Thus, it is very important to feed 
forward low-level infonnation to higher levels (bottom up) as early as possible to prevent 
wastage. 

C. FABRICATION PROCESS FLOW - BASIC STEPS [12] 

Each processing step requires certain areas be defined on the chip by the appro¬ 
priate masks during the fabrication flow. The integrated circuit is viewed as a set of pat¬ 
terned layers of doped silicon, polysilicon, metal and insulating silicon dioxide. In com¬ 
mon practice, a layer must be patterned before the next layer of material is applied on the 
chip. Lithography is the process used to transfer a pattern to a layer on the chip. Since 
each layer has its own distinct patterning requirements, by using a different mask, it is 
necessary to repeat the lithographic sequence for every layer. 

To illustrate the fabrication steps involved in patterning silicon dioxide through 
optical lithography, first examine the process flow shown in Figure 14. 
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Figure 14. Process Steps Required for Patterning of Silicon Dioxide (From [12]) 
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The thermal oxidation of the silicon surface is the process beginning, where an 
oxide layer of about one micrometer thickness is created on the substrate (Figure 14(b)). 
The entire oxide surface is then covered with a layer of photoresist that is essentially a 
light-sensitive, acid-resistant organic polymer, initially insoluble in the developing solu¬ 
tion (Figure 14(c)). The exposed areas become soluble so they are no longer resistant to 
etching solvents if the photoresist material is exposed to ultraviolet (UV) light. To expose 
the photoresist selectively, it is necessary to cover some of the areas on the surface with a 
mask during exposure. When the structure with the mask on top is exposed to UV light, 
areas covered by the opaque features on the mask are shielded. In the areas where the UV 
light can pass through, on the other hand, the photoresist is exposed and becomes soluble 
(Figure 14(d)). 

The type of photoresist that is initially insoluble and then becomes soluble after 
exposure to UV light is called positive photoresist. The process sequence shown in Figure 
14 uses positive photoresist. There is another type of photoresist that is initially soluble 
and becomes insoluble (hardened) after exposure to UV light, called negative photoresist. 
If negative photoresist is used in the photolithography process, the areas which are not 
shielded from the UV light by the opaque mask features become insoluble, whereas the 
shielded areas can subsequently be etched away by a developing solution. Negative 
photoresists are more sensitive to light, but their photolithographic resolution is not as 
high as that of the positive photoresists. Therefore, the manufacturing of high-density in¬ 
tegrated circuits use negative photoresits less commonly. 

Following the UV exposure step, a solvent can remove the unexposed portions of 
the photoresist. Therefore, the silicon dioxide regions not covered by hardened photore¬ 
sist can be etched away either by using a chemical solvent (HF acid) or by using a dry 
etch (plasma etch) process (Figure. 14(e)). Note that at the end of this step, an oxide win¬ 
dow is obtained that reaches down to the silicon surface (Figure 14(f)). It is now possible 
to strip the remaining photoresist from the silicon dioxide surface by using another sol¬ 
vent, leaving the patterned silicon dioxide feature on the surface as shown in Figure 
14(g). 
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The sequence of process steps illustrated in detail in Figure 14 actually accom¬ 
plishes a single pattern transfer onto the silicon dioxide surface, as shown in Figure 15. 
The fabrication of semiconductor devices requires several such pattern transfers to be 
performed on silicon dioxide, polysilicon, and metal. The basic patterning process used in 
all fabrication steps, however, is quite similar to the one shown in Figure 14. Also note 
that for the accurate generation of high-density patterns required in sub-micron devices, 
electron beam (E-beam) lithography is used instead of optical lithography. 



Figure 15. The Result of a Single Lithographic Patterning Sequence on Silicon Diox¬ 
ide, without Showing the Intermediate Steps. (From [12]) 

D. LAYOUT DESIGN RULES [11] 

The physical mask layout of any circuit to be manufactured using a particular 
process must confonn to the layout design rules. These rules usually specify the mini¬ 
mum allowable line widths for physical objects on-chip such as metal and polysilicon in¬ 
terconnects or diffusion areas, minimum feature dimensions, and minimum allowable 
separations between two such features. If a metal line width is made too small, for exam¬ 
ple, it is possible for the line to break during the fabrication process or afterwards, result¬ 
ing in an open circuit. If two lines are placed too close to each other in the layout, they 

may form an unwanted short circuit by merging during or after the fabrication process. 
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The main objective of design rules is to achieve a high overall yield and reliability while 
using the smallest possible silicon area for any circuit manufactured with a particular 
process. 

Note that there is usually a trade-off between higher yield, which is obtained 
through conservative geometries, and better area efficiency, obtained through aggressive, 
high-density placement of various features on the chip. The layout design rules which are 
specified for a particular fabrication process normally represent a reasonable optimum 
point in terms of yield and density. It must be emphasized, however, that the design rules 
do not represent strict boundaries which separate “correct” designs from “incorrect” ones. 
A layout which violates some of the specified design rules may still result in an opera¬ 
tional circuit with reasonable yield, whereas another layout observing all specified design 
rules may result in a circuit which is not functional and/or has very low yield. In general, 
observing the layout design rules will significantly increase the probability of fabricating 
a successful product with high yield. 

The design rules are usually described in two ways. 

• Micron rules, in which the layout constraints such as minimum feature 
sizes and minimum allowable feature separations, are stated in terms of 
absolute dimensions in micrometers, or, 

• Lambda rules, which specify the layout constraints in terms of a single pa¬ 
rameter and, thus, allow linear, proportional scaling of all geometrical 
constraints. 

E. LAYOUT OF CMOS INVERTER [12] 

The circuit consists of one nMOS and one pMOS transistor. First, it is necessary 
to create the individual transistors according to the design rules. Assume an attempt to 
design the inverter with minimum-size transistors. The minimum diffusion contact size 
then detennines the width of the active area, which is necessary for source and drain con¬ 
nections, and the minimum separation from diffusion contact to both active area edges. 

The width of the polysilicon line over the active area, which is the gate of the 
transistor, is typically taken as the minimum poly width (Figure 16). Then, the overall 
length of the active area is simply detennined by the following sum: (minimum poly 
width) + 2 (minimum poly-to- contact spacing) + 2 (minimum spacing from contact to 
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active area edge). It is necessary to place the pMOS transistor in an n-well region and the 
pMOS active area, and the minimum n-well overlap over n+ dictate the minimum size of 
the n-well. 

The distance between the nMOS and the pMOS transistor is detennined by the 
minimum separation between the n+ active area and the n-well (Figure 17). The polysili¬ 
con gates of the nMOS and the pMOS transistors are usually aligned. The final step in the 
mask layout is the local interconnections in metal, for the output node and for the VDD 
and GND contacts (Figure 18). Notice that in order to be biased properly, the n-well re¬ 
gion must also have a VDD contact. 


minimum width 
of polysilicon 



Figure 16. Design Rule Constraints which Detennine the Dimensions of a Minimum- 

Size Transistor (From [12]) 
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F. LAYOUT OF CMOS NAND AND NOR GATES [12] 

The mask layout designs of the CMOS NAND and NOR gates follow the general 
principles examined earlier for the CMOS inverter layout. Figure 19 shows the layouts of 
a two-input NOR gate and a two-input NAND gate, using a single-layer polysilicon and 
a single-layer metal. 
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Figure 19. Layouts of a CMOS NOR2 Gate and a CMOS NAND2 Gate. (From [10]) 
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Here, the p-type diffusion area for the pMOS transistors and the n-type diffusion 
area for the nMOS transistors are aligned in parallel to allow simple routing of the gate 
signals with two parallel polysilicon lines running vertically. Also notice that the two 
mask layouts show a very strong symmetry, because the NAND and the NOR gate have a 
symmetrical circuit topology. 
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IV. BICMOS COMPOSITE OPERATIONAL AMPLIFIER 
SIMULATION RESULTS 


BiCMOS technology can make op-amp improvements possible by adding the ca¬ 
pabilities of bipolar transistors to standard CMOS designs. Since bipolar transistors have 
a larger transconductance than their CMOS equivalents, they can produce op-amps with 
a larger gain bandwidth product (GBWP), arguably the most important op-amp charac¬ 
teristic. Bipolar transistors are limited in that they do not have high input impedance like 
their CMOS equivalents. In addition, BiCMOS technology requires more advanced proc¬ 
ess techniques leading to additional steps in the fabrication process. [3] 

In this chapter, the author will investigate and improve the performance of the Lee 
Configuration B op-amp [7], standard op-amp. The first part of the chapter will cover the 
basic information on the composite op-amp design. The second part of the chapter covers 
the details of the simulation of the newly designed composite amplifier, C20A2. This sets 
the basic foundation for future students to incorporate the newly designed composite am¬ 
plifier into the GIC filter to further enhance the filter performance. The following chapter 
investigates in detail the Non-Ideal performance of GIC filter, with effects from the vari¬ 
ous parameters. 

A. BACKGROUND 

Composite amplifiers were developed in 1980. Initial investigations were centered 
on increasing the GBWP and decreasing the passive and active sensitivity of a single op- 
amp. Active compensation was examined and applied to the design of active filter net¬ 
works. The resulting composite devices had three external tenninals which resembled the 
input and output terminals of a single op-amp. [16] 

In the design of two op-amp composite amplifiers (C20A), a nullator-norator 
pairing [14] was used which yielded 136 possible combinations of op-amps. Four of the 
136 possible combination yielded acceptable results based on the following criteria: 

• The non-inverting and inverting open loop gains of each of the 136 C2A’s 
were to have no change in sign in the denominator polynomial coeffi¬ 
cients, satisfying the necessary but not sufficient conditions for stability. 
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Also, the need for a single op-amps with matched GBWP’s was to be 
eliminated resulting in low component sensitivity. 

• The external three terminal performance of the C20As was to resemble, as 
close as possible, the terminal performance of the single op-amp. 

• No right half plane zeros due to the single op-amp pole were allowed in 
the closed loop gains of the C20As (to minimize the phase shift). 

• The C20A had to have minimum gain and phase deviation from the ideal 
transfer function and extended frequency operation to justify the increased 
number of op-amps. 

Based on these criteria, the following designs in Figures 20 to 23 were found to 
yield the highest response and most accurate results. 




38 





























Figure 22. C20A3 (From [15]) 



B. OPEN LOOP GAIN, 3 DB FREQUENCY AND Q P FUNCTION 

The open loop gain input output relationships for these four C20As are shown in 
Table 5. Table 6 shows the C20A 3 dB frequency and Qp functions. Finally, Table 7 
shows the stability criteria for C20As. 
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where a is the internal resistor ratio and u D is the three dB point 


Table 6. 


C20A 3dB Frequency and Qp Functions (From [15]) 



Table 7. Stability Criteria for C20As (From [15]) 
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C. COMPOSITE OPERATIONAL AMPLIFIER BANDWIDTH 

IMPROVEMENT 

The bandwidth of the single op-amp is finite. The single op-amp implemented as 
a finite gain amplifier has a bandwidth that decreases by a factor inversely proportional to 
the gain of the circuit. In contrast to the single op-amp, the bandwidth of the composite 
op-amps can be made to decrease by a factor inversely proportional to the square root of 
the gain response. This bandwidth improvement is measured for a maximally flat gain re¬ 
sponse. This maximally flat gain response is achieved when the value of Q is equal to 
0.707. [15] The theoretical frequency response of a negative finite gain composite op- 
amp compared to the frequency response of the single op-amp shows an increase in 
bandwidth of about a decade for high gain applications. 

The excellent results and the superior stability properties of the C20A1 and 
C20A2 provide the most attractive configuration for a finite gain implementation from a 
bandwidth and stability and stability point of view. [15] The C20A2 will be the configu¬ 
rations implemented in this thesis. 

D. COMPOSITE OPERATIONAL AMPLIFIERS SENSITIVITIVES 

In addition to the significant bandwidth improvements, the C20A configuration 
also offers a decreased sensitivity to active and passive components. The 3 dB frequency 
and Q are functions of the GBWPs of the single op-amps and a, the value of the ratio 
between the two resistors in the composite op-amp. 

The finite gain transfer functions for the C20As have the general form of 


l + ns 

1 +b x s + b 2 s 2 


(5.1) 


where 


and 



(5.2) 



(5.3) 
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None of the a or b coefficients are realized through differences which guarantee 
the low sensitivity of the transfer functions, the 3 dB frequencies and Q to the circuit pa¬ 
rameters. The b coefficients are always positive which is a necessary element for the 
transfer functions stability. 

E. BLOCK DIAGRAM OF AUTHOR’S STANDARD AND COMPOSITE 
DESIGN USED IN THESIS 

The standard op-amp as shown in Figure 24 has an inverting and non-inverting 
input and produces an amplified version of the input signal at the output. The value of 
the positive rail Vdd and negative rail Vss governs the maximum amplifications, pro¬ 
vided by DC power sources external to the op-amp. 



Figure 24. Block Diagram of Standard Op-Amp 


The composite op-amp shown in Figure 25 consists of two lower-order op-amps, 
i.e., Lee configuration op-amp, coupled through resistors. The C20A2 consists of two- 
coupled standard op-amps. These arrangements can produce a higher GBWP and in¬ 
crease the useful operating frequency over that which a lower-order device can achieve. 
The input, V D d and Vss requirements, and output signal are of the same overall configu¬ 
ration as the standard op-amp. 
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Figure 25. Block Diagram of Composite Op-Amp C20A2 

F. SCHEMATIC DIAGRAM OF AUTHOR’S STANDARD AND 
COMPOSITE DESIGN USED IN THESIS 

Figure 26 is the circuit layout for the Lee configuration BiCMOS versions of the 
circuit. [6] 



Figure 26. Circuit Layout of Lee Configuration 
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This op-amp design was used to produce the composite op-amps C20A2. Figure 
27 shows the symbolic diagram of a C20A2 and Figure 28 shows the schematic diagram 
of a C20A2. C20A2 is the combination of two op-amps connected in series with the for¬ 
ward op-amp having a negative feedback loop. 



Figure 27. Symbolic Diagram of C20A2 



G. SIMULATION OF AUTHOR’S STANDARD AND COMPOSITE DESIGN 
USED IN THESIS 

This section simulates and compares the open loop gain, closed loop gain, slew 
rate and common rejection ratio of the composite op-amp to the standard op-amp. 
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1. Open Loop Gain 

Each op-amp was designed in Silvaco EXPERT then extracted and tested with 
Silvaco SMARTSPICE. Figures 29 and 30 show the open loop transfer function for Lee’s 
Configuration and C20A2, respectively. Lee’s configuration displays an open loop gain 
of 45dB while C20A2 registers an open loop gain of 49dB. The gain and a of the com¬ 
posite amplifier is 100 and 6, respectively. 


*Normal OP AMP Transfer Function - Open Loop 



Figure 29. Open Loop Transfer Function of Lee’s Configuration B 
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*c2a01 OP AMP Transfer Function - Open Loop 



Figure 30. Open Loop Transfer Function of C20A2 

Figure 31 shows the comparison between them. It shows that C20A2 has a higher 
open gain loop as compared to Lee’s basic configuration. All codes for the simulation 
appear in Appendix A. 


— + v(vout2) 

— 0 acl.v(vout) 



Figure 31. Open Loop Transfer Function - Comparison 
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2. Closed Loop Gain 

Both operational amplifiers were arranged in the configuration shown in Figure 
32 to analyze the closed loop gain. The gain of the amplifier is configured to 100. 


R2 

■VWV 


F.l 

WvV 


’ I® krfi 
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Figure 32. Closed Loop Gain Configuration 


Figure 33 shows the closed loop gain of 100 for Lee’s configuration. It has a 3dB 
comer frequency of 110 kHz. Which yielded an equivalent GBWP of about 11 MHz. 

On the other hand, Figure 34 shows the C20A2 3dB comer frequency of 1.2 
MHz, for the same gain, with a 22.4 MHz OdB crossing frequency. The corresponding 
GBWP has improved to to the equivalent of 120 MHz. All codes for the simulation ap¬ 
pear in Appendix B. 


48 



Gain (dB) 


LeeOpamp Closed Loop Gain 



Frequency (Hz) 


Figure 33. Closed Loop Gain of Lee’s Configuration 

C20A2 Closed Loop Gain 



Figure 34. Closed Loop Gain of C20A2 

3. Slew Rate 


A pulse was injected into the amplifier to observe the output response. Figure 35 
shows the slew rate of Lee’s configuration. It exhibits a slew rate of 20MV/sec. 
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Figure 35. Slew Rate of Lee’s Configuration 

Figure 36 shows the slew rate of 40MV/sec for C20A2. The slew rate has in¬ 
creased by 20MV/sec in C20A2. All codes for the simulation appear in Appendix C. 


• • tmu 1 
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Figure 36. Slew Rate of C20A2 


4 Common Mode Rejection Ratio (CMRR) 

CMRR is the ability of an operational amplifier to reject the common-mode sig¬ 
nal. It is expressed in terms of dB. Figure 37 shows a comparison between Lee’s configu- 
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ration and C20A2. It shows that C20A2 has a higher CMRR. All codes for the simulation 
appear in Appendix D. 





Figure 37. CMRR of Lee’s Configuration and C20A2 

H. SUMMARY OF SIMULATION PERFORMANCE INDICTOR 

Table 8 summarizes the performance indictor for the composite op-amp. The ad¬ 
dition of a second op-amp increases the number of poles at lower frequencies and drives 
down the curve of the transfer function at a steeper rate than Lee’s configuration op-amp. 
This steeper curve allows the frequency of the 3dB point to be pushed higher, but yet still 
reduces the gain below OdB before the impact of parasitic capacitances occur at higher 
frequencies. If these parasitic capacitances come into effect before the gain drops below 
OdB, the op-amp may saturate due to positive feedback created by phase shifts and any 
small amount of noise. It is very clear that the C20A2 is a better op-amp because of im¬ 
proved performance and was easy to implement. 



Open Loop 
Gain (dB) 

3dB Corner 
Frequency 
(Hz) ' 

CMRR 

(dB) 

Slew Rate 
(V/sec) 

GBWP (Hz) 

Lee’s configu¬ 
ration 

45 

110k 

79 

20 M 

10M 

C20A2 
(new design) 

49.42 

1.2M 

93 

40M 

12M 


Table 8. Summary of Performance Indictor 
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V. GIC FILTER SIMULATION RESULTS 


This chapter investigates a continuous analog generalizedd impedance converter 
(GIC) 4 th -order filter band pass filter in detail. Varies parameter like network sensitivity, 
effects of resistor and capacitor values and GBWP on center frequency (fo ) and Q factor 
(Qp ) are investigated. P-SPICE and MATLAB software are used to develop and investi¬ 
gate the simulation differences from the theory result for both the Ideal and Non-Ideal 
GIC Band Pass filter. Due to the complexity of the 4 th order transfer function, MAPLE 8 
software is used to derive the Transfer Function for the Ideal and Non-Ideal Op-amp for 
the GIC filter. All the results from theory, MATLAB and P-SPICE are used to compare 
the actual circuit buildup prior to future GIC filter chip fabrication. 

As mentioned in Chapter II, the GIC filter is well known for its stability and low 
sensitivity. Each GIC biquad employs two op-amps and has good high frequency per¬ 
formance. The previous chapter proved that the composite op-amp is better than the 
standard op-amplifier. The performance of the GIC filter will then be further improved 
by utilizing the newly designed composite op-amp that will not be covered in this thesis. 
A. DESIGN DETAILS 

The fourth-order filter is achieved by cascade realization. Two matching second- 
order filter blocks are connected in series (cascaded) to construct a fourth-order band 
pass filter. Although alternate design methods such as direct synthesis or follow-the- 
leader-feedback (FLF) may potentially yield slightly better perfonnance, the simplicity 
and ease of troubleshooting make cascade realization the optimal method for designing 
this filter. The Band Pass stages are gain adjustable. All GIC stages have equal capacitor 
values, unless a capacitor is required to adjust the gain. 

Figure 38 (schematic) and Figure 39 (circuit layout) depict the circuit topology. 
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Figure 39. Circuit Layout 
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B. IDEAL AND NON IDEAL TRANSFER FUNCTIONS 

The transfer function for the ideal and non-ideal op-amp for GIC have been de¬ 
rived using Maple 8 software. The initial parameters are set as center frequency (fo ) of 
50KHz and Q factor (Qp ) of 1. Refer to Appendix E for the derivation of the ideal trans¬ 
fer function for the 4th-order GIC band pass filter. 


4th-order TF (ideal): 


T(s) = 


{QpRC) • G + (2Q p RC)s + (2 Q p 2 + l)s 2 + 


k RCj 


5 + 


v RC 2 j 


(5.1) 


For the purposes of investigating the GIC band pass filter behaviors, resistor and 
capacitor values are chosen. Substituting design values (R= 1 kQ, Qp = 10, C = 3.183 
nF) produce the transfer function below: 


T(s) = 


_ (9.86e8)s 2 _ 

7 + (6.28£>4> 3 + (1.98<?1 1)7 + (6.20el 5)5 + 9.74£>21' 


(5.2) 


4th-order TF (non-ideal) = (with design values substituted into T(s)): 

^ , (9.87e22)/+(2.04e30)/+(1.05e37)i 2 

1 (a) -- 

/ +(2.13e7)F +(2.27el4)/ +(1.21e21)i 5 +(3.28e27).f +(4.06e32)i 3 +(5.67e38).y 2 +(2.62e42).v+2.43e48 

(5.3) 

C. SENSITIVITIES ANALYSIS FOR S T ( S ) R , S T( s) C , S ro o R AND S QP R 

Since the 4th-order filter is achieved by cascade realization, for the ease of calcu¬ 
lation, the sensitivities are calculated based on the 1 st stage (2 nd order) and multiple a fac¬ 
tor of 2 for the 2 nd stage. 

1. S T(s , r = S(Numerator, R) - S(Denominator, R) 

= R( ( d (Num)/ d R)/Num -(dd (Deno)/ d R)/Deno (5.4) 
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where 


T(s)= 2 


/ s 1 \ 

rcQp s 2 + -^ + ^-7 

^ rcQp r 2 c 2 J 


(5.5) 


8 ( s 

8 (Num)/ 8 r = | 2 


= -2 


dr y r c Qp 
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r cQp 


(5-6) 

(5.7) 


8 (Deno)/ 8 r = — 
dr 
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r c Qp ~ 2 ~ 2 


r~ c 


Qp(s 2 r 2 c 2 -l) 

~~ s 2 r 2 c 2 Qp + s r c + Qp 


(5-8) 

(5-9) 


(5.10) 

(5-11) 


Therefore, S T ( S ) K for the 2 nd stage 


-2 


Qp(s 2 r 2 c 2 -l) 
s 2 r 2 c 2 Qp + s r c + Qp 


(5.12) 


2. S T( s) = S(Num,C) - S(Deno,C) 

= C( (8 (Num)/ 8 C)/Num -(8 (Deno)/ G C)/Dcno) 


8 (_ s 

8 (Num)/ 3 c = ^r\ 2 


8c \ r c Qp 


-2 


rc Qp 
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(5.13) 

(5-14) 

(5.15) 


























(5.16) 


where 


8 (Deno)/ 8 c = 


St( S) C - 
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8c 
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r c Qp + r 2 c 2 


(5.17) 


(5.18) 


Qp{s 2 r 2 c 2 - 1) 
s 2 r 2 c 2 Qp 


(5.19) 


Therefore, St (S ) C for the 2 nd stage 

= 2 Qp(s 2 r 2 c 2 - 1) 

s 2 r 2 c 2 Qp + s r c + Qp 


3. S(oo R = S(Num,r) - S(Deno,r) 

= r((5 (Num)/ 6 r)/Num -(8 (Deno)/ 6 r)/Dcno) 


co 0 = 1 /RC 
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= 0 
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(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 


Therefore, S^o 11 for the 2 nd stage 

= -2 
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4. S Qp r = S(Num,R) - S(Deno,R) 

= R( (8 (Num)/ 8 R)/Num - ( 8 (Deno)/ 8 R)/Deno) 

(5.27) 


where 


Qp = Rq/R 

8 rq 

8 (Num)/ 5 R = ^ - 


(5.28) 

(5.29) 
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(5.30) 

(5.31) 
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Therefore, Sq p r for the 2 nd stage 


(5.32) 

(5.33) 


(5.34) 


PLOTTING THE IDEAL AND NON-IDEAL TRANSFER FUNCTION 
USING MATLAB AND P-SPICE SOFWARE 


Please refer to Appendix F for the ideal GIC filter MATLAB program and Ap¬ 
pendix G for the non-ideal GIC filter MATLAB program. 


For the ideal case, the fo , Gain and Qp are set as 50KHz, 12 dB and Qp = 10, re¬ 
spectively. Figure 40 shows the MATLAB Ideal transfer function plot at the second stage 
output. Figure 41 shows the MATLAB non-ideal transfer function plot at the second 
stage output. 
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Figure 40. MATLAB Ideal Transfer Function 


Second Stage Non-Ideal Band Pass - Closeup scale fcx Mag and Phase 
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Figure 41. MATLAB Non-Ideal Transfer Function 


From the MATLAB simulation, the results are as follows: 
Non-Ideal (1 st stage): fo =46.951 kHz, Gain=10.9 dB, Qp =9.516 
Ideal (1 st stage): fo =49.815 kHz, Gain=12 dB, <2p=T0.097 
Non-Ideal (2 nd stage): fo =46.951 kHz, Gain=10.9 dB, Qp=17.353 
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Ideal (2 nd stage): fo =49.975 kHz, Gain=12 dB, Qp=\5.1 

Figure 42 shows the P-SPICE Frequency response plot. From the P-SPICE simu¬ 
lation, the results are as follows: 

fo = 45.394 kHz, Gain= 11.997 dB, Qp =16.091 


P Sch BP UK 



Figure 42. P-SPICE Frequency Response Plot 


1. Comparing the Ideal Values with MATLAB Program Values 

With the ideal simulation, the fo is closer to 50 kFIz in the 2 nd stage as compared 

to the 1 st stage. Flowever, the Qp has increased from 10.097 to 15.7 in the 2 nd stage. 
When compared to the non-ideal case, the fo is slightly lower than those in the ideal 
case second stage. However, the Qp is 10% higher than in the ideal case 2 nd stage. Both 
in MATFAB are quite close. 

2. Comparing the MATLAB Program Values and P-SPICE Values 

With non-ideal simulation, the fo has decreased to 45.395 kHz. The Qp has also 
dropped slightly to 16.01. The main reason is due mainly to the operation of this fire- 
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quency near the GBWP of the op-amp. The P-SPICE shows the actual values to expect 
in real circuit construction. 

*Note: Reference is taken from the non-ideal MATLAB program values or all the 
remaining comparisons. 

E. ANALYSIS OF THE EFFECT OF 10% VARIATION ON VALUE OF 

RESISTOR ON fo AND Qp OF THE NON-IDEAL TRANSFER 

FUNCTION 

Please refer to Appendix H for the non-ideal MATLAB program (non-ideal case 
with 10% increase in resistor value). The MATLAB simulation for the non-ideal case 
with no resistor value and capacitor value changes is as follows: 

Non-Ideal: fo =46.951 kHz, Gain= 10.9, <2p = 17.353 

Ligure 43 shows the non-ideal MATLAB plot (non-ideal case with 10% increase 
in resistor value). The MATLAB simulation (non-ideal case with 10% increase in resis¬ 
tor value) is as follows. 

Non-Ideal: fo =46.951 kHz, Gain= 10.9, Qp=17.353 


Second Stage Non-Ideal Band Pass with R 10% Increase • Closeup scale for Mag and Phase 



Ligure 43. MATLAB Simulation (Non-Ideal Case with 10% Increase in Resistor 
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Value) 


Figure 44 shows the P-SPICE Frequency response plot. From the P-SPICE simu¬ 
lation, the results are as follows. 

fo =41.687 kHz, Gain=l 1.877, Qp =14.499 



Figure 44. P-SPICE Frequency Response Plot (with 10% Increase in Resistor Value) 

1. Comparing the Values in Both MATLAB Programs (with and without 
10% Increase in Resistor Value) 

With 10% increase in Resistor value, there is no difference between the gain, and 
there is no difference in fo and Qp . The results show that an increase of 10% in Resistor 
value does not affect Gain, Qp and fo for this circuit. 

2. Comparing the MATLAB Simulation Program Values (with 10% In¬ 
crease in Resistor Value) and P-SPICE Values (with 10% Increase in 
Resistor Value) 

With a 10% increase in resistor value, the fo has dropped to 41.687 kHz, Qp has 
also dropped to 14.499. However, the gain has increased to 11.877. 
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F. ANALYSIS OF THE EFFECT OF 10% VARIATION ON VALUE OF 

CAPACITOR ON fo AND Qp OF THE NON-IDEAL TRANSFER 

FUNCTION 

Please refer to Appendix I for the non-ideal MATLAB program (non-ideal case 
with 10% increase in capacitor value). The MATLAB simulation for the non-ideal case 
with no resistor value and capacitor values changes is as follows. 

Non-Ideal: fo =46.951 kHz, Gain= 10.9, <2p=17.353 

Figure 45 shows the non-ideal MATLAB plot (non-ideal case with 10% increase 
in capacitor value). The MATLAB simulation (non-ideal case with 10% increase in ca¬ 
pacitor value) is as follows. 

Non-Ideal: fo =42.971 kHz, Gain=l 1.1, Qp=16.875 


Second Stage Nan Ideas Band Pass wtfti C 10% rarest* - Ooseup scale tor Mag and Phase 



Figure 45. MATLAB Simulation (Non-Ideal Case with 10% Increase in Capacitor 

Value) 

Figure 46 shows the P-SPICE Frequency response plot. From the P-SPICE simu¬ 
lation, the results are as follows. 

fo =41,687 kHz, Gain=12.053, Qp=16.595 
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Figure 46. P-SPICE Frequency Response Plot (with 10% Increase in Capacitor 

Value) 

1. Comparing the Values in Both MATLAB Programs (with and without 
10% Increase in Capacitor Value) 

With 10% increase in capacitor value, the gain has increased slightly to 11.1, both 
fo and Qp have decreased to 42.97kHz and 16.875, respectively. In summary, with an 
increase of capacitor value, fo and Qp will decrease and gain will increase. 


2. Comparing the MATLAB Program Values (with 10% Increase in 

Capacitor Value) and P-SPICE Values (with 10% Increase in Capaci¬ 
tor Value) 

With 10% increase in capacitor value, the fo has dropped to 41.687kHz, Qp has 
also dropped slightly to 16.595. However, the gain has increased to 12.053. 

G. ANALYSIS OF THE EFFECT OF 10% VALUE INCREASE IN RESISTOR 
AND 50% DECREASE IN GBWP OF THE NON-IDEAL TRANSFER 
FUNCTION 

Please refer to Appendix J for the non-ideal MATLAB program (non-ideal case 
with 10% increase in resistor value and 50% decrease of GBWP). The MATLAB simula- 
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tion for the non-ideal case with a 10% increase in resistor value and 100% of GBWP is 
as follows. 

Non-ideal case with 10% increase in resistor value, 100% GBWP: fo =46.951 
kHz, Gain=10.9, gp=17.353 

Refer to Figure 47 for the MATLAB plot (non-ideal case with 10% increase in 
resistor value, 50% GBWP). The MATLAB simulation for the non-ideal case with 10% 
increase in resistor value with 50% decrease of GBWP is as follows. 

Non-Ideal case: fo =44.563 kHz, Gain=9.18, <2p=15.56 


Socond Stag* Non kJoal Band Paaa wtlh R 10% Inoraaaa ft S0% QOWP Cloaaup aoala lor Mag and Pfcaao 



Figure 47. MATLAB Plot (Non-Ideal Case with 10% Increase in Resistor Value and 

50% Decrease of GBWP) 


Figure 48 shows the P-SPICE Frequency response plot. From the P-SPICE simu¬ 
lation, the results are as follows. 


fo =41.687 kHz, Gain=l 1.877, Qp =14.83 
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Figure 48. P-SPICE Frequency Response Plot (with 10% Increase in Resistor Value 

and 50% Decrease of GBWP) 

1. Comparing the Values in Both MATLAB Programs (with 10% In¬ 
crease in Resistor Value for 100% and 50% Reduction of GBWP) 

With 50% GBWP reduction, the gain has reduced to 9.18, the fo is reduced to 
44.563 kHz, and Qp is reduced to 15.56. In summary, with 50% GBWP reduction, it will 
reduce the gain, fo and Qp. 

2. Comparing the MATLAB Program Values (with 10% Increase in Re¬ 
sistor Value for 50% Reduction of GBWP) and P-SPICE Values 
(with 10% Increase of Resistor Value for 50% Reduction of GBWP) 

With 50% GBWP reduction, the fo and Qp has reduced slightly while the gain 
has increased slightly to 11.877. 

H. ANALYSIS OF THE EFFECT OF 10% VALUE INCREASE IN 

CAPACITOR AND 50% DECREASE IN GBWP OF THE NON-IDEAL 
TRANSFER FUNCTION 

Please refer to Appendix K for the non-ideal MATLAB program (non-ideal case 
with 10% increase in capacitor value and 50% decrease of GBWP). The MATLAB simu¬ 
lation for the non-ideal case with 10% increase in capacitor value and 100% of GBWP 
are as follows. 

Non-Ideal case: fo =46.951 kHz, Gain=10.9, Qp =17.353 
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Refer to Figure 49 for the MATLAB plot (non-ideal case with 10% increase in 
capacitor value, 50% GBWP). The MATLAB simulation for the non-ideal case with a 
10% increase in capacitor value and 50% of GBWP are as follows. 

non-ideal case with 10% increase in capacitor value, 50% GBWP: 
fo =40.90ckHz, Gain=9.53, <2p = 16.06 


Srcond Stag* Non IdMl Hand llil with C 10% Inaaaaa and 50% o< QDWP Clean.* acala to Mag and PTiaaa 



Figure 49. MATLAB Plot (Non-Ideal Case with 10% Increase in Capacitor Value 

and 50% Decrease of GBWP). 

Figure 50 shows the P-SPICE frequency response plot. From the P-SPICE simu¬ 
lation, the results are as follows. 

fo =41.687 kHz, Gain=12.053, Qp= 16.522 
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Figure 50. P-SPICE Frequency Response Plot (with 10% Increase in Capacitor 

Value and 50% Decrease of GBWP) 

1. Comparing the Values in Both MATLAB Programs (with 10% In¬ 
crease in Capacitor Value for 100% and 50% Reduction of GBWP) 

With 50% GBWP reduction, the gain has reduced to 9.53, the fo has reduced to 
40.90kHz, and Qp has reduced to 16.06. In summary, with 50% GBWP reduction, it will 
reduce the gain, fo and Qp . 

2. Comparing the MATLAB Program Values (with 10% Increase in 
Capacitor Value for 50% Reduction of GBWP) and P-SPICE Values 
(with 10% Increase in Capacitor Value for 50% Reduction of GBWP) 

With 50% GBWP reduction, the fo , gain and Qp have reduced slightly as com¬ 
pared to P-SPICE. 

I. SUMMARY OF SIMULATION RESULTS 

There are some differences when comparing the MATLAB program ideal with 
non-ideal cases. When the operating frequency is close to ideal, the percentage errors be¬ 
come zero, i.e., toward the ideal case. 
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A 10% increase in resistor values does not affect gain, Qp and fo . However, 
with 10% increase of capacitor values, fo and Qp decrease with an increase in gain. 
With 50% GBWP reduction (in the case of 10% increase in resistor value and capacitor 
value), it reduces the gain, fo and Qp . All values from the P-SPICE program when 
compared with MATLAB program values are well within 10% error. 

After adjusting for non-ideal conditions, the final theoretical value of the circuit 
built in the laboratory performed very close to the theoretical value. The final value of the 
center frequency was 50.560 kHz compared to the target value of 50.0 kHz (1.1% 
higher). The value of Qp was 9.711 compared to a target of 10.0, which is less than 3% 
error. 

J. PRACTICAL CIRCUIT ANALYSIS 

In the process of practical circuit building, it was observed that none of the mod¬ 
els, i.e., MATLAB and P-SPICE, exactly matches the actual, but P-SPICE shows closer 
values. However, replacing the Rq and C to achieve the said fo of 50 kHz did several it¬ 
eration processes. Please refer to Table 9 for the results. 



MATLAB Ideal 

MATLAB non-ideal 

Rq 

C 

fo 

Qp 

fo 

Qp 

10k 

3.183n 

49.975k 

15.7 

46.951k 

17.353 

6.667k 

3.183n 

49.656k 

10.4 

46.792k 

10.138 

6.667k 

2.86n 

53.386k 

10.54 

51.725k 

10.156 

6.667k 

3.01 n 

52.521k 

10.645 

49.338k 

10.33 



P-SPICE 

Actual 

Rq 

C 

fo 

Qp 

fo 

Qp 

10k 

3.183n 

45.395k 

16.091 

47.9 

14.5 

6.667k 

3.183n 

45.394k 

10.584 

47.9 

9.39 

6.667k 

2.86n 

50.119k 

10.605 

52.5 

8.89 

6.667k 

3.01 n 

47.863k 

10.655 

50.560k 

9.711 


Table 9. Summary of Component Values 


All the practical values used in the circuit are also calculated as shown in Table 


10 . 
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3 100 

3U150 2651 

1 00E-07 

31 33 

316.31 

213.27 

hiu 

3141592654 

1.006‘01 

JlfcSi 

3163. f0 

2132.68 

O.Olu 

3H1E912B64 

2.006-06 

mis 

Ifi4l .5-5 

'056 3+ 

0.02u 

3U1S0.2ES4 

1.506-07 

21.22 

21 £21 

limits 

0.15U 

31+159 255+ 

2,20€-07 

14.47 

1+4.190 

95.94 

Q.Km 

3H1S4SSS4 

S.B0&0T 

12.73 

197 3? 

65.31 

0.2fili 

31+153.2654 

7.B06-10 

4244.13 

42441.22 

26425 58 

750 a 

31+159.2654 

6.5JE-10 

+581 03 

46810.28 

31952.49 

6HHp 

3141&9.2S&3 

G.B0E-10 

80B4.il 

sso+ios 

3505350 

E80O 

314169.25&i 

+.70E-10 

6772.55 

67725.51 

4537409 

47CO 

3141=$,R5+ 

3.3PE-CS 

964.55 

9615.75 

6462.65 

2.2n 


Table 10. Summary of Calculated Values 


MATLAB, P-SPICE and actual frequency response for the practical circuit are 
depicted in the Figures 51, 52 and 53, respectively. 



Frequency (rad/sec) 


Figure 51. MATFAB Plot: Center Frequency = 49.338 kHz, Qp = 10.33 
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Figure 52. 


Figure 53. 


147.S63K,21.6921 
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P-SPICE Frequency Response: Center Frequency = 47.853 kFIz, Qp = 

10.655 


REF 11.0 dBm MARKER 50 960. O Hz 

2 dS/OlV RANGE 20.0 dBm 7.02 dBm 



_ / 

START 40 OOO. O Hz STOP BO OOO. O Hz 

RBW 300 Hz VBW 1 KHz ST . 6 SEC 


Actual Frequency Response: Center Frequency = 50.560 kFIz, Qp = 9.711 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The primary goal of this thesis was to investigate the continuous analog general¬ 
ized impedance converter (GIC) 4th-order band pass fdter, ideal and non-ideal cases, in 
detail with software aids such as MATLAB, P-SPICE and MAPLE. MATLAB software 
depicted more accurate data when compared to theoretical results whereas P-SPICE 
demonstrated the practice aspect of the circuit design. In the actual circuit setup, there are 
too many components variation like tolerates in resistor value and capacitor value. How¬ 
ever, all this can be overcome by implementing the design in VLSI. 

The gain, fo and Q factor are not affected by 10% increase in resistor values. 
However, 10% increase in capacitor value causes center frequency and Q factor to de¬ 
crease and gain to increase. With 50% GBWP reduction (for 10% increase in Resistor 
value and capacitor values), it will reduce the gain, center frequency and Q factor. All the 
P-SPICE program values when compared to MATLAB program values are all well 
within 10% error. The final design of the circuit built in the laboratory performed very 
close to the stated specifications. The final value of the center frequency was 50.560 kHz 
compared to the target value of 50.0 kHz (1.1% higher). The value of Q factor was 9.711 
compared to a target of 10.0, which is less than 3% error. 

The secondary objective was to investigate the previous student’s designed 
BiCMOS standard operational amplifier and improve it as a composite operational ampli¬ 
fier. The composite amplifier was introduced during the Analog VLSI course taught by 
Professor Michael at the Naval Postgradute School, and the author was very impressed by 
the design and has proven its worth in this thesis. The results show that the gain band¬ 
width product (GBWP), common mode rejection ratio (CMRR), and open loop gain im¬ 
proved considerably. 

The author tried to simulate the filter design initially using PSPICE and found it 
very difficult as the actual component technical data are different from the design re¬ 
quirements. Efforts were made to contact http://www.orcadpcb.com for the details. How¬ 
ever, those parameters were unable to be alttered to suit the filter design. The author sub- 
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sequently designed the composite amplifier (C20A2) in Silvaco EXPERT software, and 
then extracted and simulated it with the Silvaco SmartSpice software. It was possible to 
alter all parameters by changing the width and length of the gates using the Silvaco 
EXPERT software. 

With the improvements from the composite operational amplifier, it sets the basic 
foundation for future students to incorporate the newly designed composite operational 
amplifier into the GIC filter to enhance the filter performance further. 
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APPENDIX A. OPEN LOOP 


A. OPEN LOOP - LEE’S CONFIGURATION [6], SPICE NETLIST 

*Normal OP AMP Transfer Function - Open Loop 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ GND DC 0.0 AC 10m 

VI- VIn- GND DC 0.0 

*OPAMP CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 VOUT VDD CMOSP1 W= 150U L= 5.2U 

M6 VOUT VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl VOUT VN3 1.032P 

ckill vout gnd 1.966p 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*R1 RlinVin- lk 
*R2 VOutVin- 100k 

* Input Stimulus, Options and Measurments 
*.TRAN lp 16n On 

.AC DEC 20 10 1E+12 
.plot v(VOut) v(VDD) v(VSS) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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B. OPEN LOOP - C20A2 CONFIGURATION SPICE NETLIST 

*c2a01 OP AMP Transfer Function - Open Loop 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ GND DC 0.0 AC 10m 

VI- Vin- GND DC 0.0 

*OPAMP_l CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 vin+2 VDD CMOSP1 W= 150U L= 5.2U 

M6 Vin+2 VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl vin+2 VN3 1.032P 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*OPAMP_2 CIRCUIT 

Q12 VN12 VIN-2 VN52 NPN1 AREA=2.560P 
Q22 VN22 VIN+2 VN52 NPN1 AREA=2.560P 
M32 VN52 VN42 VSS VSS CMOSN1 W= 27.2U L= 5.2U 
M42 VN42 VN42 VSS VSS CMOSN1 W= 13.2U L= 5.2U 
M52 VDD GND VN42 vss CMOSN1 W= 5.2U L= 33.2U 
M92 VDD VN2 VOUT2 VDD CMOSP1 W= 150U L= 5.2U 
M62 VOUT2 VN42 VSS VSS CMOSN1 W= 67.2U L= 5.2U 
C12 VOUT2 VN32 1.032P 

Ml02 VN32 VSS VN22 VDD CMOSP1 W= 5.2U L= 5.2U 
M72 VN12 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 
M82 VN22 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 

* interconnect bet opl and op2, alpha=6.10 lk 10k 
Rlc vin-2 gnd lk 

R2c vin-2 vout2 8k 
*feedback R 
*R1 rlinVin- lk 
*R2 VOut2 Vin- 100k 
ckill vin+2 gnd 1.966p 
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ckill2 vout2 gnd 1.966p 

*Input Stimulus, Options and Measurments 

*.TRAN lp 16n On 

.AC DEC 20 10 1E+12 

.plot v(VOut2) 

.plot v(VDD) 

.plot v(VSS) 
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APPENDIX B. CLOSED LOOP 


A. CLOSED LOOP - LEE’S CONFIGURATION [6], SPICE NETLIST 

*Normal OP AMP Transfer Function - closed loop 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ GND DC 0.0 AC 10m 

VI- R1 in GND DC 0.0 

*OPAMP CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 VOUT VDD CMOSP1 W= 150U L= 5.2U 

M6 VOUT VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl VOUT VN3 1.032P 

ckill vout gnd 1.966p 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
R1 RlinVin- lk 
R2 VOut Vin- 100k 

* Input Stimulus, Options and Measurments 
*.TRAN lp 16n On 

.AC DEC 20 10 1E+12 
.plot v(VOut) v(VDD) v(VSS) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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B. CLOSED LOOP - C20A2 CONFIGURATION SPICE NETLIST 

*c2a01 OP AMP Transfer Function -closed Loop 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ GND DC 0.0 AC 10m 

VI-R1 in GND DC 0.0 

*OPAMP_l CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 vin+2 VDD CMOSP1 W= 150U L= 5.2U 

M6 Vin+2 VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl vin+2 VN3 1.032P 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*OPAMP_2 CIRCUIT 

Q12 VN12 VIN-2 VN52 NPN1 AREA=2.560P 
Q22 VN22 VIN+2 VN52 NPN1 AREA=2.560P 
M32 VN52 VN42 VSS VSS CMOSN1 W= 27.2U L= 5.2U 
M42 VN42 VN42 VSS VSS CMOSN1 W= 13.2U L= 5.2U 
M52 VDD GND VN42 vss CMOSN1 W= 5.2U L= 33.2U 
M92 VDD VN2 VOUT2 VDD CMOSP1 W= 150U L= 5.2U 
M62 VOUT2 VN42 VSS VSS CMOSN1 W= 67.2U L= 5.2U 
C12 VOUT2 VN32 1.032P 

Ml02 VN32 VSS VN22 VDD CMOSP1 W= 5.2U L= 5.2u 
M72 VN12 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 
M82 VN22 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 

* interconnect bet opl and op2, alpha=6.10 lk 10k 
Rlc vin-2 gnd lk 

R2c vin-2 vout2 8k 
*feedback R 
R1 rlinVin- lk 
R2 VOut2 Vin- 100k 


80 



ckill vin+2 gnd 1.966p 

ckill2 vout2 gnd 1.966p 

*Input Stimulus, Options and Measurments 

*.TRAN lp 16n On 

.AC DEC 20 10 1E+12 

.plot v(VOut2) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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APPDENDIX C. SLEW RATE 


A. SLEW RATE - LEE’S CONFIGURATION [6], SPICE NETLIST 

^Normal OP AMP Transfer Function - slew rate 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ vin+ gnd PULSE(0 5 0 0.0 IpS 0.0 lps lOmS ) 

VI- VIn- vout 0.0 

*OPAMP CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 VOUT VDD CMOSP1 W= 150U L= 5.2U 

M6 VOUT VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl VOUT VN3 1.032P 

ckill vout gnd 1.966p 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*R1 RlinVin- lk 
*R2 VOutVin- 100k 

* Input Stimulus, Options and Measurments 
*.TRAN lm 20m On 

.TRAN 10PS 400NS INS 
*.AC DEC 20 10 1E+6 
.plot v(VOut) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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B. SLEW RATE - C20A2 CONFIGURATION SPICE NETLIST 

*c2a01 OP AMP Transfer Function -slew rate 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ vin+ gnd PULSE(0 5 0 0.0 IpS 0.0 lps lOmS ) 

VI- VIn- vout2 0.0 

*OPAMP_l CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 vin+2 VDD CMOSP1 W= 150U L= 5.2U 

M6 Vin+2 VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl vin+2 VN3 1.032P 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*OPAMP_2 CIRCUIT 

Q12 VN12 VIN-2 VN52 NPN1 AREA=2.560P 
Q22 VN22 VIN+2 VN52 NPN1 AREA=2.560P 
M32 VN52 VN42 VSS VSS CMOSN1 W= 27.2U L= 5.2U 
M42 VN42 VN42 VSS VSS CMOSN1 W= 13.2U L= 5.2U 
M52 VDD GND VN42 vss CMOSN1 W= 5.2U L= 33.2U 
M92 VDD VN2 VOUT2 VDD CMOSP1 W= 150U L= 5.2U 
M62 VOUT2 VN42 VSS VSS CMOSN1 W= 67.2U L= 5.2U 
C12 VOUT2 VN32 1.032P 

Ml02 VN32 VSS VN22 VDD CMOSP1 W= 5.2U L= 5.2u 
M72 VN12 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 
M82 VN22 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 

* interconnect bet opl and op2, alpha=6.10 lk 10k 
Rlc vin-2 gnd lk 

R2c vin-2 vout2 8k 
*feedback R 
*R1 rlinVin- lk 
*R2 VOut2 Vin- 100k 
ckill vin+2 gnd 1.966p 
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ckill2 vout2 gnd 1.966p 

*Input Stimulus, Options and Measurments 

*.TRAN lm 20m On 

.TRAN 10PS 400NS INS 

*.AC DEC 20 10 1E+6 

.plot v(VOut2) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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APPENIDX D. CMRR 


A. CMRR - LEE’S CONFIGURATION [6], SPICE NETLIST 

^Normal OP AMP Transfer Function - CMRR 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ Vin- DC 0.0 AC lOOOp 

VI- Vin- GND DC 0.0 

*VI+ VIn+ GND DC 0.0 AC lOOOp 

♦VI- Vin- GND DC 0.0 AC lOOOp 

*OPAMP CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 VOUT VDD CMOSP1 W= 150U L= 5.2U 

M6 VOUT VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl VOUT VN3 1.032P 

ckill vout gnd 1.966p 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*R1 Rlin Vin- lk 
*R2 VOutVin- 100k 

* Input Stimulus, Options and Measurments 
*.TRAN lp 16n On 

.AC DEC 20 10 1E+6 
.plot v(VOut) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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B. CMRR - C20A2 CONFIGURATION SPICE NETLIST 

*c2a0 I OP AMP Transfer Function -CMRR 

.INCLUDE AMIparams.txt 

*Power Supplies 

VDD VDD GND DC 5.0 

VSS VSSGND DC-5.0 

* Signal 

VI+ VIn+ vin- DC 0.0 ac lOOOp 

VI- Vin- GND DC 0.0 

*OPAMP_l CIRCUIT 

.MODEL CMOSP1 PMOS 

.MODEL CMOSP2 PMOS 

.MODEL CMOSN1 NMOS 

.MODEL CMOSN2 NMOS 

.MODEL NPN1 NPN 

Q1 VN1 VIN-VN5 NPN1 AREA=2.560P 

Q2 VN2 VIN+ VN5 NPN1 AREA=2.560P 

M3 VN5 VN4 VSS VSS CMOSN1 W= 27.2U L= 5.2U 

M4 VN4 VN4 VSS VSS CMOSN1 W= 13.2U L= 5.2U 

M5 VDD GND VN4 vss CMOSN1 W= 5.2U L= 33.2U 

M9 VDD VN2 vin+2 VDD CMOSP1 W= 150U L= 5.2U 

M6 Vin+2 VN4 VSS VSS CMOSN1 W= 67.2U L= 5.2U 

Cl vin+2 VN3 1.032P 

M10 VN3 VSS VN2 VDD CMOSP1 W= 5.2U L= 5.2U 
M7 VN1 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
M8 VN2 VN1 VDD VDD CMOSP1 W= 30U L= 5.2U 
*OPAMP_2 CIRCUIT 

Q12 VN12 VIN-2 VN52 NPN1 AREA=2.560P 
Q22 VN22 VIN+2 VN52 NPN1 AREA=2.560P 
M32 VN52 VN42 VSS VSS CMOSN1 W= 27.2U L= 5.2U 
M42 VN42 VN42 VSS VSS CMOSN1 W= 13.2U L= 5.2U 
M52 VDD GND VN42 vss CMOSN1 W= 5.2U L= 33.2U 
M92 VDD VN2 VOUT2 VDD CMOSP1 W= 150U L= 5.2U 
M62 VOUT2 VN42 VSS VSS CMOSN1 W= 67.2U L= 5.2U 
C12 VOUT2 VN32 1.032P 

Ml02 VN32 VSS VN22 VDD CMOSP1 W= 5.2U L= 5.2u 
M72 VN12 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 
M82 VN22 VN12 VDD VDD CMOSP1 W= 30U L= 5.2U 

* interconnect bet opl and op2, alpha=6.10 lk 10k 
Rlc vin-2 gnd lk 

R2c vin-2 vout2 8k 
*feedback R 
*R1 rlinVin- lk 
*R2 VOut2 Vin- 100k 
ckill vin+2 gnd 1.966p 
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ckill2 vout2 gnd 1.966p 

*Input Stimulus, Options and Measurments 

*.TRAN lp 16n On 

.AC DEC 20 10 1E+6 

.plot v(VOut2) 

.plot v(VDD) 

.plot v(VSS) 

.end 
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APPENDIX E. DERIVATION OF IDEAL TRANSFER FUNCTION 
FOR 4TH-ORDER GIC BANDPASS FILTER 
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APPENDIX F. IDEAL GIC FILTER MATLAB PROGRAM 


%Chcong,Hcng Wan 

Syms s 

R=1E3; 

G=l/R; 

Qp=5; 

f=10 A 4; 

wp=2*pi*f; 

C=(R*wp) A -l; 

s=tf('s'); 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

T1=(Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8)/(Y1*Y4*(Y5+Y6)+Y2*Y3*(Y7+Y8)); 

T2=(Y1*Y4*Y5+Y1*Y5*Y8+Y2*Y3*Y7- 

Y1*Y6*Y7)/(Y1*Y4*(Y5+Y6)+Y2*Y3*(Y7+Y8)); 

T3=(Y1*Y4*Y5+Y2*Y3*Y7)/(Y1*Y4*(Y5+Y6)+Y2*Y3*(Y7+Y8)); 

T1 

[num 1 ,den 1 ]=tfdata(T 1, V) 

figure(l) 

step(Tl) 

figure(2) 

bode(Tl) 
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APPENDIX G. NON- IDEAL GIC FILTER MATLAB PROGRAM 


%Cheong,Heng Wan 
%Non-ideal GIC opamp 
clear all; 

Syms s 

%parameters 

R=le03; 

G=l/R; 

Qp=i; 

f=20e03; 

wp=2*pi*f 

C=(R*wp) A -l; 

s=tf('s'); 

wtl=le07 %as given 
wt2=le07 %as given 
Al=wtl/s; 

A2=wt2/s; 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

Tl= ((Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8) + 

((Y7*(Y1+Y3)*(Y2+Y5+Y6))/A1)) / ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + 

(Y1 *(Y 4+Y7+Y 8)*(Y2+Y 5+Y 6))/A 1 + (Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

T2= ((Y2*Y3*Y7+Y1*Y5*(Y4+Y8)-Y1*Y6*Y7) + (Y5*(Y1+Y3)*(Y4+Y7+Y8))/A2 ) 
/ ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + (Y1*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A1 + 
(Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

[num 1 ,den 1 ]=tfdata(T 1 ,'v'); 

%Full scale for Mag and phase 
figure(l); 

Bode(Tl); 

title('Non Ideal Band Pass - Full scale for Mag and Phase'); 
grid; 

%Closeup scale for Mag and phase 
figure(2); 

Bode(Tl); 

title('Non Ideal Band Pass - Closeup scale for Mag and Phase'); 
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axis([l.le05,1.4e05,0,6.2]); 

grid; 

%Closeup scale for Mag only 
[MAG 1,PHASE 1, W1 ] = BODE(Tl); 

MAGDB1= 20*logl0(MAGl); 
figure(3); 

semilogx(W 1 ,MAGDB 1 

title('Non Ideal Band Pass - Closeup scale for Mag only’); 

xlabel(’W in log scale’); 

ylabel('Mag in DB’); 

axis([l.le05,1.4e05,0,6.2]); 

grid; 

MAGDBmax=max(M AGDB1) 

[Y,I] = max(MAGDBl); 

Wlmax = W1(I,1) 
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APPENDIX H. NON-IDEAL MATLAB PROGRAM (NON IDEAL 
CASE WITH 10% INCREASE IN R) 


%effect of 10% R on Wo,Qp 
%Cheong,Heng Wan 
%Non-ideal GIC opamp 
clear all; 
clc; 

Syms s 

%parameters 

R=(le03)*l.l; 

G=l/R; 

Qp=10; 

f=20e03; 

wp=2*pi*f 

C=(R*wp) A -l; 

s=tf('s'); 

wtl=le07 %as given 
wt2=le07 %as given 
Al=wtl/s; 

A2=wt2/s; 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

Tl= ((Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8) + 

((Y7*(Y1+Y3)*(Y2+Y5+Y6))/A1)) / ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + 

(Y1 *(Y 4+Y7+Y 8)*(Y2+Y 5+Y 6))/A 1 + (Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

T2= ((Y2*Y3 *Y7+Y1 * Y5 *(Y4+Y8)-Y 1 *Y6*Y7) + (Y5*(Y1+Y3)*(Y4+Y7+Y8))/A2 ) 
/ ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + (Y1*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A1 + 
(Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

[num 1 ,den 1 ]=tfdata(T 1,'v'); 

%Closeup scale for Mag and phase 
figure(l); 

Bode(Tl); 

title('Non Ideal Band Pass'); 
axis([l.le05,1.4e05,0,6.2]); 
grid; 
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%Closeup scale for Mag only 
[MAG 1,PHASE 1, W1 ] = BODE(Tl); 

MAGDB1= 20*logl0(MAGl); 

MAGDBmax=max(MAGDB 1) 

[Y,I] = max(MAGDBl); 

Wlmax = W1(I,1) 

MAFDB3 db=MAGDBmax-3 
%Non-ideal GIG opamp - reference 
%parameters 
Ra=le03; 

Ga=l/Ra; 

Qpa=10; 
fa=20e03; 
wpa=2*pi*fa 
Ca=(Ra*wpa) A -l; 
s=tf('s'); 

wtla=le07 %as given 
wt2a=le07 %as given 
Ala=wtla/s; 

A2a=wt2a/s; 

Yla=Ga; 

Y2a=Ga; 

Y3a=s*Ca; 

Y4a=Ga; 

Y5a=0; 

Y6a=Ga; 

Y7a=Ga/Qpa; 

Y8a=s*Ca; 

Tla= ((Yla*Y4a*Y5a+Y3a*Y7a*(Y2a+Y6a)-Y3a*Y5a*Y8a) + 
((Y7a*(Yla+Y3a)*(Y2a+Y5a+Y6a))/Ala)) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 
(Yla*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/Ala + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
T2a= ((Y2a*Y3a*Y7a+Yla*Y5a*(Y4a+Y8a)-Yla*Y6a*Y7a) + 
(Y5 a* (Y1 a+Y3 a) * (Y4a+Y7a+Y8a))/A2a ) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 

(Y1 a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A 1 a + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
[num 1 a, den 1 a]=tfdata(T 1 a, V); 

%Closeup scale for Mag and phase 
figure(2); 

Bode(Tla,Tl); 

title('Non Ideal Band Pass - 10% variable in R'); 
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axis([l.le05,1.4e05,0,6.2]); 

grid; 

%Closeup scale for Mag only 

[MAG 1 a,PHASE 1 a, W1 a] = BODE(Tla); 

MAGDBla= 20*logl0(MAGla); 

MAGDBmaxa=max(MAGDB 1 a) 

[Ya,Ia] = max(MAGDBla); 

W1 maxa = W1 a(Ia, 1) 

MAFDB3 dba=M AGDBmaxa-3 
%Differences 

Diffindb = abs(MAGDBmaxa-MAGDBmax) 

Diff_in_db_percent = abs(((MAGDBmaxa-MAGDBmax)/MAGDBmaxa)* 100 ) 
Diff_in_w = abs(Wlmaxa-Wlmax) 

Diff_in_w_percent = abs(((W 1 maxa-W 1 max)/W 1 maxa)* 100) 
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APPENDIX I. NON-IDEAL MATLAB PROGRAM (NON IDEAL 
CASE WITH 10% INCREASE IN C) 


%effect of 10% C on Wo,Qp 
%Cheong,Heng Wan 
%Non-ideal GIC opamp 
clear all; 
clc; 

Syms s 

%parameters 

R=(le03); 

G=l/R; 

Qp=10; 
f=20e03; 
wp=2*pi*f 
C=((R*wp) A -l)* 1.1; 
s=tf('s'); 

wtl=le07 %as given 
wt2=le07 %as given 
Al=wtl/s; 

A2=wt2/s; 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

Tl= ((Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8) + 

((Y7*(Y1+Y3)*(Y2+Y5+Y6))/A1)) / ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + 

(Y1 *(Y 4+Y7+Y 8)*(Y2+Y 5+Y 6))/A 1 + (Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

T2= ((Y2*Y3 *Y7+Y1 * Y5 *(Y4+Y8)-Y 1 *Y6*Y7) + (Y5*(Y1+Y3)*(Y4+Y7+Y8))/A2 ) 
/ ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + (Y1*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A1 + 
(Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

[num 1 ,den 1 ]=tfdata(T 1 ,'v'); 

%Closeup scale for Mag and phase 
figure(l); 

Bode(Tl); hold on; 
title('Non Ideal Band Pass'); 
axis([0.95e05,1.4e05,-l,6.2]); 
grid; 
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%Closeup scale for Mag only 
[MAG 1,PHASE 1, W1 ] = BODE(Tl); 

MAGDB1= 20*logl0(MAGl); 

MAGDBmax=max(MAGDB 1) 

[Y,I] = max(MAGDBl); 

Wlmax = W1(I,1) 

MAFDB3 db=MAGDBmax-3 
%Non-ideal GIG opamp - reference 
%parameters 
Ra=le03; 

Ga=l/Ra; 

Qpa=10; 
fa=20e03; 
wpa=2*pi*fa 
Ca=(Ra*wpa) A -l; 
s=tf('s'); 

wtla=le07 %as given 
wt2a=le07 %as given 
Ala=wtla/s; 

A2a=wt2a/s; 

Yla=Ga; 

Y2a=Ga; 

Y3a=s*Ca; 

Y4a=Ga; 

Y5a=0; 

Y6a=Ga; 

Y7a=Ga/Qpa; 

Y8a=s*Ca; 

Tla= ((Yla*Y4a*Y5a+Y3a*Y7a*(Y2a+Y6a)-Y3a*Y5a*Y8a) + 
((Y7a*(Yla+Y3a)*(Y2a+Y5a+Y6a))/Ala)) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 
(Yla*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/Ala + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
T2a= ((Y2a*Y3a*Y7a+Yla*Y5a*(Y4a+Y8a)-Yla*Y6a*Y7a) + 
(Y5 a* (Y1 a+Y3 a) * (Y4a+Y7a+Y8a))/A2a ) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 

(Y1 a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A 1 a + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
[num 1 a, den 1 a]=tfdata(T 1 a, V); 

%Closeup scale for Mag and phase 
figure(2); 

Bode(Tla,Tl); 

title('Non Ideal Band Pass - 10% variable in C'); 
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axis([0.95e05,1.4e05,-l ,6.2]); 
grid; 

%Closeup scale for Mag only 

[MAG 1 a,PHASE 1 a, W1 a] = BODE(Tla); 

MAGDBla= 20*logl0(MAGla); 

MAGDBmaxa=max(MAGDB 1 a) 

[Ya,Ia] = max(MAGDBla); 

W1 maxa = W1 a(Ia, 1) 

MAFDB3 dba=M AGDBmaxa-3 
%Differences 

Diffindb = abs(MAGDBmaxa-MAGDBmax) 

Diff_in_db_percent = ab s (( (M AGD B maxa—MA GD B max)/MAGD B maxa) * 100 ) 
Diff_in_w = abs(Wlmaxa-Wlmax) 

Diff_in_w_percent = abs(((W 1 maxa-W 1 max)/W 1 maxa)* 100) 
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APPENDIX J. NON-IDEAL MATLAB PROGRAM (NON-IDEAL 
CASE WITH 10% INCREASE IN R AND 50% DECREASE OF 

GBWP). 


% effect of 10% R on Wo,Qp - 50% of GBWP 

%Chcong,Hcng Wan 

%Non-ideal GIC opamp 

clear all; 

clc; 

Syms s 

%parameters 

R=(le03)*l.l; 

G=l/R; 

Qp=10; 

f=20e03; 

wp=2*pi*f 

C=(R*wp) A -l; 

s=tf('s'); 

wtl=(le07)*0.5 %as given 
wt2=(le07)*0.5 %as given 
Al=wtl/s; 

A2=wt2/s; 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

11= ((Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8) + 

((Y7*(Y1+Y3)*(Y2+Y5+Y6))/A1)) / ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + 

(Y1 *(Y 4+Y7+Y 8)*(Y2+Y 5+Y6))/A 1 + (Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

12= ((Y2*Y3*Y7+Y1*Y5*(Y4+Y8)-Y1*Y6*Y7) + (Y5*(Y1+Y3)*(Y4+Y7+Y8))/A2 ) 
/ ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + (Y1*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A1 + 

(Y3 * (Y 4+Y 7+Y 8) * (Y2+Y 5+Y 6 ))/A2 + 

((Y1 + Y 3 ) * (Y 4+Y 7+Y 8) * (Y 2+Y 5+Y 6))/( A1 * A2)); 

[num 1 ,den 1 ]=tfdata(T 1,'v'); 

%Closeup scale for Mag and phase 
figure(l); 

Bode(Tl); 

title('Non Ideal Band Pass'); 
axis([l.le05,1.4e05,0,6.2]); 
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grid; 

%Closeup scale for Mag only 
[MAG 1,PHASE 1, W1 ] = BODE(Tl); 

MAGDB1= 20*logl0(MAGl); 

MAGDBmax=max(M AGDB1) 

[Y,I] = max(MAGDBl); 

Wlmax = W1(I,1) 

MAFDB3 db=MAGDBmax-3 
%Non-ideal GIG opamp - reference 
%parameters 
Ra=le03; 

Ga=l/Ra; 

Qpa=10; 
fa=20e03; 
wpa=2*pi*fa 
Ca=(Ra*wpa) A -l; 
s=tf('s'); 

wtla=(le07) %as given 
wt2a=(le07) %as given 
Ala=wtla/s; 

A2a=wt2a/s; 

Yla=Ga; 

Y2a=Ga; 

Y3a=s*Ca; 

Y4a=Ga; 

Y5a=0; 

Y6a=Ga; 

Y7a=Ga/Qpa; 

Y8a=s*Ca; 

Tla= ((Yla*Y4a*Y5a+Y3a*Y7a*(Y2a+Y6a)-Y3a*Y5a*Y8a) + 
((Y7a*(Yla+Y3a)*(Y2a+Y5a+Y6a))/Ala)) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 

(Y1 a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A 1 a + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
T2a= ((Y2a*Y3a*Y7a+Yla*Y5a*(Y4a+Y8a)-Yla*Y6a*Y7a) + 
(Y5a*(Y1 a+Y3a)*(Y4a+Y7a+Y8a))/A2a ) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 
(Yla*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/Ala + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
[num 1 a, den 1 a]=tfdata(T 1 a, V); 

%Closeup scale for Mag and phase 
figure(2); 

Bode(Tla,Tl); 
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title('Non Ideal Band Pass - 10% variable in R & 50% of GBWP’); 

axis([l.le05,1.4e05,0,6.2]); 

grid; 

%Closeup scale for Mag only 

[MAG 1 a,PHASE 1 a, W1 a] = BODE(Tla); 

MAGDBla= 20*logl0(MAGla); 

MAGDBmaxa=max(MAGDB 1 a) 

[Ya,Ia] = max(MAGDBla); 

Wlmaxa = Wla(Ia,l) 

MAFDB3 dba=MAGDBmaxa-3 
%Differences 

Diffindb = abs(MAGDBmaxa-MAGDBmax) 

Diff_in_db_percent = ab s (( (M AGDB maxa—MA GD B max)/MAGD Bmaxa) * 100 ) 
Diff_in_w = abs(Wlmaxa-Wlmax) 

Diff_in_w_percent = abs(((W 1 maxa-W 1 max)/W 1 maxa)* 100) 
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APPENDIX K. NON-IDEAL MATLAB PROGRAM (NON-IDEAL 
CASE WITH 10% INCREASE IN C AND 50% DECREASE OF 

GBWP). 


% effect of 10% C on Wo,Qp & 50% of GBWP 

%Chcong,Hcng Wan 

%Non-ideal GIC opamp 

clear all; 

clc; 

Syms s 

%parameters 

R=(le03); 

G=l/R; 

Qp=10; 
f=20e03; 
wp=2*pi*f 
C=((R*wp) A -l )* 1.1; 
s=tf('s'); 

wtl=(le07)*0.5 %as given 
wt2=(le07)*0.5 %as given 
Al=wtl/s; 

A2=wt2/s; 

Y1=G; 

Y2=G; 

Y3=s*C; 

Y4=G; 

Y5=0; 

Y6=G; 

Y7=G/Qp; 

Y8=s*C; 

11= ((Y1*Y4*Y5+Y3*Y7*(Y2+Y6)-Y3*Y5*Y8) + 

((Y7*(Y1+Y3)*(Y2+Y5+Y6))/A1)) / ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + 

(Y1 *(Y 4+Y7+Y 8)*(Y2+Y 5+Y6))/A 1 + (Y3*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A2 + 
((Y1+Y3)*(Y4+Y7+Y8)*(Y2+Y5+Y6))/(A1*A2)); 

12= ((Y2*Y3*Y7+Y1*Y5*(Y4+Y8)-Y1*Y6*Y7) + (Y5*(Y1+Y3)*(Y4+Y7+Y8))/A2 ) 
/ ((Y2*Y3*(Y7+Y8)+Y1*Y4*(Y5+Y6)) + (Y1*(Y4+Y7+Y8)*(Y2+Y5+Y6))/A1 + 

(Y3 * (Y 4+Y 7+Y 8) * (Y2+Y 5+Y 6 ))/A2 + 

((Y1 + Y 3 ) * (Y 4+Y 7+Y 8) * (Y 2+Y 5+Y 6))/( A1 * A2)); 

[num 1 ,den 1 ]=tfdata(T 1,'v'); 

%Closeup scale for Mag and phase 
figure(l); 

Bode(Tl); hold on; 
title('Non Ideal Band Pass'); 
axis([0.95e05,1.4e05,-l,6.2]); 
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grid; 

%Closeup scale for Mag only 
[MAG 1,PHASE 1, W1 ] = BODE(Tl); 

MAGDB1= 20*logl0(MAGl); 

MAGDBmax=max(M AGDB1) 

[Y,I] = max(MAGDBl); 

Wlmax = W1(I,1) 

MAFDB3 db=MAGDBmax-3 
%Non-ideal GIG opamp - reference 
%parameters 
Ra=le03; 

Ga=l/Ra; 

Qpa=10; 
fa=20e03; 
wpa=2*pi*fa 
Ca=(Ra*wpa) A -l; 
s=tf('s'); 

wt 1 a= 1 e07 %as given 
wt2a=le07 %as given 
Ala=wtla/s; 

A2a=wt2a/s; 

Yla=Ga; 

Y2a=Ga; 

Y3a=s*Ca; 

Y4a=Ga; 

Y5a=0; 

Y6a=Ga; 

Y7a=Ga/Qpa; 

Y8a=s*Ca; 

Tla= ((Yla*Y4a*Y5a+Y3a*Y7a*(Y2a+Y6a)-Y3a*Y5a*Y8a) + 
((Y7a*(Yla+Y3a)*(Y2a+Y5a+Y6a))/Ala)) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 

(Y1 a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A 1 a + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
T2a= ((Y2a*Y3a*Y7a+Yla*Y5a*(Y4a+Y8a)-Yla*Y6a*Y7a) + 
(Y5a*(Y1 a+Y3a)*(Y4a+Y7a+Y8a))/A2a ) / ( 
(Y2a*Y3a*(Y7a+Y8a)+Yla*Y4a*(Y5a+Y6a)) + 
(Yla*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/Ala + 
(Y3a*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/A2a + 
((Yla+Y3a)*(Y4a+Y7a+Y8a)*(Y2a+Y5a+Y6a))/(Ala*A2a)); 
[num 1 a, den 1 a]=tfdata(T 1 a, V); 

%Closeup scale for Mag and phase 
figure(2); 

Bode(Tla,Tl); 
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title('Non Ideal Band Pass - 10% variable in C & 50% of GBWP'); 

axis([0.95e05,1.4e05,-l ,6.2]); 

grid; 

%Closeup scale for Mag only 

[MAG 1 a,PHASE 1 a, W1 a] = BODE(Tla); 

MAGDBla= 20*logl0(MAGla); 

MAGDBmaxa=max(MAGDB 1 a) 

[Ya,Ia] = max(MAGDBla); 

Wlmaxa = Wla(Ia,l) 

MAFDB3 dba=MAGDBmaxa-3 
%Differences 

Diffindb = abs(MAGDBmaxa-MAGDBmax) 

Diff_in_db_percent = abs(((MAGDBmaxa-MAGDBmax)/MAGDBmaxa)*100 ) 
Diff_in_w = abs(Wlmaxa-Wlmax) 

Diff_in_w_percent = abs(((W 1 maxa-W 1 max)/W 1 maxa)* 100) 
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